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We present a 3D photoionization model of the planetary nebula NGC 6302, one of 

the most complex and enigmatic objects of its kind. It's highly bipolar geometry and 

dense massive disk, coupled with the very wide range of ions present, from neutral 

^^ I species up to Si*^ , makes it one of the ultimate challenges to nebular photoionization 

modelling. 

Our MOCASSIN model is composed of an extremely dense geometrically thin cir- 
cumstellar disk and a large pair of diffuse bipolar lobes, a combination which was 
O ■ necessary to reproduce the observed emission-line spectrum. The masses of these com- 

^ I ponents, 2.2 M© and 2.5 Mq respectively, gives a total nebular mass of 4.7 M©, of 

C/5 , which 1.8 Mq (39%) is ionized. Discrepancies between our model fit and the observa- 

tions are attributed to complex density inhomogeneities in the nebula. The potential 
to resolve such discrepancies with more complex models is confirmed by exploring a 
range of models introducing small-scale structures. Compared to solar abundances he- 
^ ' lium is enhanced by 50%, carbon is slightly subsolar, oxygen is solar, and nitrogen is 

^^ I enhanced by a factor of 6. These all imply a significant 3rd dredge-up coupled with hot- 

'/^ ■ bottom burning CN-cycle conversion of dredged-up carbon to nitrogen. Aluminium is 

*^ ' also depleted by a factor of 100, consistent with depletion by dust grains. 

The central star of NGC 6302 is partly obscured by the opaque circumstellar disk, 
[*--. ' which is seen almost edge-on, and as such its properties are not well constrained. 

^^ I However, emission from a number of high-ionization 'coronal' lines provides a strong 

constraint on the form of the high-energy ionizing flux. We model emission from the 
central star using a series of stellar model atmospheres, the properties of which are 
constrained from fits to the high-ionization nebular emission lines. Using a solar abun- 
dance stellar atmosphere we are unable to fit all of the observed line fiuxes, but a 
^% ' substantially better fit was obtained using a 220,000 K hydrogen-deficient stellar at- 

H I mosphere with log g — 7.0 and L^, — 14, 300 L0. The H-dcficient nature of the central 

star atmosphere suggests that it has undergone some sort of late thermal pulse, and 
fits to evolutionary tracks imply a central star mass of 0.73-0.82 Mq. Timescales for 
these evolutionary tracks suggest the object left the top of the asymptotic giant branch 
~2100 years ago, in good agreement with studies of the recent mass-loss event that 
formed one pair of the bipolar lobes. Based on the modelled nebular mass and central 
star mass we estimate the initial mass of the central star to be 5.5 Mq, in approximate 
agreement with that derived from evolutionary tracks. 
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1 INTRODUCTION evolutionary phase the star ejects a large fraction of its mass 

during brief phases (10'^ — 10* yrs) of high mass-loss (up to 

Planetary nebulae (PNe) are one of the last evolutionary m-t at -i\ rp,- i ' ^ n iu 

J y ' ■'10 Mq yr ). ihis mass-loss eventually exposes the cen- 

stages of the majority of low and intermediate mass stars (1- ^^^j ^^^^ ^^ ^^^ ^^^^ ^^^^ ^^.^^^^ ^^^^^^1^ .^ ^^^ Uy ^^^ 

8 Mq). During the previous asymptotic giant branch (AGB) 
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Figure 1. The planetary nebula NGC 6302 as imaged with the 
Wide Field Camera on the Hubble Space Telescope (HST). Blue is 
He II, brown is H«, and white is [S ll]. The image is approximately 
2 X 2.5 arcmin or 0.6 X 0.75 pc at the distance of NGC 6302. Image 
copyright: NASA and ESA. 



ionizes the surrounding nebula making it visible for study. 
Despite supposed spherically symmetric mass-loss while on 
the AGB, many PNe show point- or axi-symmetric shapes 
and complex morphologies that have been interpr eted as evi- 
dence for the influence of a binary companion (e.g. IPe Marcol 
120091 ) ■ This and other outstanding issues in PNe research 
such as their small-scale s tructures (e.g. iGoncalves et al.l 
I2OO1I : iMatsuura et al.ll2009l ) and discrepancies between for- 
bidden li ne and recombination line abundances for heavy 
eleme nts (|Liu et al.ll200(]| : iTsamis et al.|[20o3 : 1 Wesson et all 
I2OO5I ) suggest that out understanding of the late stages of 
stellar evolution is far from complete. Despite this, PNe are 
still regularly used as standard candles, metallicity indica- 
tors and tracers of stell ar populations in distant galaxies 
(e.g. iBuzzoni et al.ll2006l ). If PNe are to be used in this way 
our understanding of their structure and evolution must be 
greatly improved, both through studies of large samples o f 
such objects (e.g. IViironen et al.l l2009l: Sabin et al. 201CI), 
their precursors the AGB stars (e.g. van d er Veen fc Habind 
[l988; Zijlstra et al. 2001; Wright et al. 200i), and through 
individual c ase studies to understand the most complex ex- 
amples fe.g. iGongalves et aLJlJOOel : IWareing et al.ll2006l ). 

NGC 6302 is one of the most complex and extreme ex- 
amples of a planetary nebula. The nebula is broadly bipo- 
lar, with a second narrower pair of lobes visible further out 
(see Figure [l| . It has a highl y pinched waist that is typi- 
cal of 'butterfly' bipolar PNe JBalick fc Frankll2002l ). as is a 
highly complex structure with many clumps and knots. The 
bipolar structure has been attributed to the presence of a 
very dense circumstellar disk that is seen almost edge-on 
and almost completely obscures the centra l star {Av ~ 6.6, 
iMatsuura et aklboosl : ISzvszka et al.l bOOgl ) . which has only 
recently been detected. The disk radia tes strongly in the 
infrared (e.g. iLester fc DincrsteinI 11984 ) with evidence for 



significant fractions of very large grains (JHoare et al.lll992l ) 



and crystalline silicates JMolster et al.ll200ll '). NGC 6302 has 
b een classified as a T ype I PN, both according to the criteria 
or lPeimbert fc Torres-Pciinbert (1983, N/0 ^ 0.5 and H e/H 
> 0.125) and the criteria of lKingsburgh fc Barlowl (|l994l ') for 
a Milky- Way Type I PN, that N/O ^ 0.8. The high helium 
and nitrogen abundances observed in the nebula imply a po- 
tentially massive progenitor, while observ ations of emission 
lines from very high ionization stages (e.g. Ashlev fc Hvlandl 
1 19881 : [Rowlands et al.|[T993 : ICasassus et al.|[200ol ') suggest an 
extremely hot central star. 

Attempts to understand NGC 6302 have been hindered 
by its complexity and the inability to study its central star 
and th erefore determine its evolutionary history. lAUer et al.l 
(|l981f) could find "no uniform density theoretical model that 
could reproduce the line intensiti es in NGC 6302 in any sat- 
isfactory fashion", while Lame fc Ferlandl (|l991<) could not 
produce a physically plausible mo del for the nebula wi thout 
requiring shock-excit ation, which lOriglia et al.l (|1993l ) and 
ICasassus et al.l (|200d ') found strong evidence against. How- 
ever, many of these models were forced to assume spherical 
symmetry due to the computational constraints of full 3- 
dimensional modeling, yet the inhomogeneous morphology 
of NGC 6302 clearly requires a 3-dimensional model if reli- 
able conclusions are to be drawn. 

In this paper we present a 3D photoionization model 
of NGC 6302 using the 3D Monte Carl o photoioniza- 
tion and radiativ e transfer code MOCASSIN (jErcolano et al.l 
|2003| , I2OO5I . |2008| ). The model presented here is a fully 3- 
dimensional model of NGC 6302 that includes both gas and 
dust, though we will only discuss the gas photoionization 
model in this paper. We will present the dust radiative trans- 
fer model in a future paper. In Section 2 we discuss the 
changes made to the MOCASSIN code to enable the modeling 
of this object and outline the observations that our models 
are to reproduce. In Section 3 we describe the model used 
and the parameters used to build the model. In Section 4 
we describe the modeling process and how different quan- 
tities varied according to changes in difi^erent parameters. 
In Section 5 we present the model results, including fits to 
the observed emission-line spectrum and the ionization and 
temperature structure of the nebula. Finally, in Section 6 
we describe the resulting nebula properties and discuss the 
evolutionary history of the object. 



2 THE MOCASSIN CODE 

MOCASSIN is a 3D photoionization and radiative transfer 
code that uses Monte Carlo met hods to solve the rad iative 
transfer in a photoionized nebula IjErcolano et al.ll2003l ) . The 
Monte Carlo approach allows it to fully treat all physical 
processes in 3D, self-consistently treating the diffuse com- 
ponent of the radiation field that ID photoionization codes 
are unable to deal with. The code is highly versatile and can 
model any structure or morphology, introducing unlimited 
subgrids to resolve detail in areas of the model. It also al- 
lows multiple non-central ionizing sources and multiple gas 
chemistries. MOCASSIN also has a complete radiative trans- 
fer dust model including the full treatment of all dust pro- 
cesses such as photoe lectric heating and gas-grain collisions 
(jErcolano et alj 20051). The atomic data base included opac- 
ity data from IVerner et al.l l|l993i) and IVerner fc YakovlevI 
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l|l995h . energy levels, collisions strengths and transition 
prob a .bilities from Version 5.2 of the CHIANtQ (|Dere et all 
1 19971 : iLandi et al.ll2006i . and references therein) and the im- 
proved hj^drogen_aBd2ieliunifi;ee-bound continuous emission 
data of lErcolano fc Storevl(|2006r ). In this section we describe 
a number of changes made to the code as part of this work, 
as well as details of the computational facilities that the 
models were run on, and the observations that were to be 
reproduced by the models. 

2.1 2-dimensional simulations in MOCASSIN 

MOCASSIN is a 3-dimensional code, but many objects ex- 
hibit cylindrical symmetry along a rotation axis. Such ob- 
jects may effectively be modeled in two dimensions, r and 
z, greatly reducing the computational and memory require- 
ments of large and complicated models. We have devel- 
oped and tested a 2-dimensional version of the MOCASSIN 
code, which integrates the observations over 6 from to 360 
(where 6 is the angle between two directions in the z = 
plane). Since the nebula is symmetrical upon reflection in 
the z = Q plane, only positive z values need be modeled, 
reducing the processing time further. This approximation 
still properly treats the diffuse component of the radiation 
field and thorough benchmarking has been performed that 
ful ly reproduces a ll of th e mocassin benchmarks described 
in lErcolano et al.l l|2003r ). Although the models described 
here were run in 2D, we will often discuss the resulting 3D 
structure of the nebula, since the code reproduces the full 
3D structure in the outputs. The final model was run in 3D 
to confirm the validity of the 2D modeling approach and no 
differences were found. The new r outines are included in the 
latest public version of the code (|Ercolano et aLlbOOST ). 

2.2 Computational details 

The models presented in this work were run on the 
hiperspace facilities at University College London, includ- 
ing the Enigma SUN Microsystems V880 multiprocessor 
computer, the Keter SUN Microsystems V880 and V890 mi- 
croprocessor computers and the Dell Legion Cluster, which 
uses Intel's dual-core technology. Simulations were typically 
run with 8-16 processors and approximately 32 GB of mem- 
ory. The necessary computational processing time varied de- 
pending on the complexity of the model and the required 
level of convergence of the Monte Carlo simulation, typically 
an iteration-to-iteration variation of the fraction of neutral 
hydrogen of <1% in >60% of the cells in the model. 

2.3 Observations to reproduce 

MOCASSIN simulations produce several outputs that are used 
to compare the model results with observations. The outputs 
can also be adapted to the exact setup of the observations 
by simulating the slit or aperture through which the obser- 
vations were made and the inclina tion angle of the nebula 
to the line of sight. As shown by IConcalves et al.l (|2006l ) 
this is particularly important for long-slit observations of 
highly inhomogeneous nebulae such as NGC 7009 and the 

-*- http://www.ukssdc.ac.uk/solar/chianti/ 
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broad range of different measurements for NGC 6302 of the 
same emission line is testament to this. We have searched 
the literature for observations of NGC 6302 that span all 
wavelength regimes and that include lines from all the im- 
portant elements and their ionization stages. The list in Ta- 
ble [T] was chosen based on spectral coverage and resolution. 
Where measurements overlap the most precise measurement 
was taken, or that which was accompanied by complemen- 
tary information such as electron temperatures or densities 
derived from diagnostic line ratios. 

The UV spectrum is mainly taken from the Interna- 
tional Ultrav i olet E xplorer (lUE) spectrum presented by 
iTsamis et al.l l|2003r ) and c omplemented by observations of 
the [Ne v] 3426 A Une by iRowlands et alj l|l994 ) for com- 
parison with the [Ne v] lines in the mi d-IR. The ma i n opt i- 
cal emission-line spectrum was that of iTsamis et al.l (|2003), 
a deep optical spectrum made with a fixed slit aligned 
along the polar axis of the nebula that covers all aspects 
of the nebula. These observations include many of the di- 
agnostic line ratios that will be important to diagnose the 
density and ionization structure within the nebula. This 
main line list is supplemented by lin es from the far- r ed re - 
gion using the pu blished spectra of iDanziger et al.l (119731 ) 
and I Groves et al.l ([2002). The infrared line list is split be- 
tween the Infrared Space Observatory (ISO) observations 
usi ng the Short Wavelength Spe ctrometer (SWS) presented 
by iBeintema fc PottaschI 1119991) and t he Lo ng Wavelength 
Spectrometer (LWS) from iLiu et al.l (120011). and comple- 
mented by deep echelle spectroscopy from ICasassus et al.l 
(2000). These include many fine-structure lines in the mid- 
and far-infrared which have the advantage of only having a 
weak temperature dependence due to their low excitation 
energies. 

In cases where the observed lines are known to be blends 
that are inseparable in the observations (e.g. the blend of 
[Ar iv] A4711 and [Ne iv] A4712), we accounted for the pres- 
ence of the undesired component by calculating its intensity 
relative to other lines of the same species using the equib 
code. 

Where possible all observations have been converted 
onto a scale relative to H/3 for comparison with mod- 
eled line fluxes. For observations where this is not pos- 
sible we compare observed fluxes, calculating model line 
fluxes us ing the distance to N GC 6302 of 1.17 kpc mea- 
sured by iMeaburn et al.l (|2008r ) . We also compare the II/3 
flux for the entire nebula from our models with that 
measured by ITsamis et a\\ (|2003t ) as F(H/3) = 2.82 x 
10~^^ erg cm~^ s~^. Using their derived reddening for 
NGC 6302 (c(H/3) = 1.44) we obtain a dereddened H/3 flux 
of 7.76x10"^° ergscm-^s"^ or L(/3) = 1.27 xlO^'^ ergs s-\ 
which can be compared with model outputs for the nebula. 



2.4 Observational constraints on the central star 
luminosity 

The central star of NGC 6302 was identifled for the first 
time in deep Hubb le Space Telescope images reported by 
ISzvszka et al.l (|2009l ). who were able to perform photometric 
measurements of the central star in two filters and there- 
fore make estimates of the surface temperature, extinc- 
tion and stellar luminosity. However, at such high temper- 
atures and extinctions, the exact solution is not well con- 
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Table 1. Spectroscopic observations to be matched by the MOCASSIN photoionisation model of NGC 6302. 



Region 


Wavelength range 


Telescope 


Lines 


Aperture / slit 


Reference 


Ultraviolet 


1150-1975 A 


lUE SWP 


7 


10.3" X 23" (oval) 


Tsamis et al. (2003 1 


Ultraviolet 


3426 A 


Mt. Lemmon 1.5-m 


1 


48" (circular) 


Rowlands et al. C1994) 


Optical 


3040-7400 A 


ESO 1.52-m 


30 


210" X 2" 


Tsamis et al. (2003) 


Optical-IR 


3300-8600 A 


Siding Springs 2.3-m 


1 


2" (long-slit) 


Groves et al. (20021 


Optical-IR 


3500-20,000 A 


Cerro Tololo 92-cm 


1 


17" X 34" 


Danziger et al. (1973) 


Near-infrared 


2.7-4.8 /im 


UKIRT (echelle) 


12 


3"x 3" 


Casassus et al. (2000) 


Mid-infrared 


2.4-45 fim 


ISO SWS 


14 


14" X 20" 


Beinterna & Pottasch (19991 


Mid-infrared 


8-13 ^m 


UKIRT N-band 


5 


4" (circular) 


Casassus et al. (2000) 


Far-infrared 


43-198 Aim 


ISO LWS 


7 


40" (circular) 


Liu et al. (20011 



strained and the authors were forced to use constraints on 
the age of the nebula, combined with evolutionary mod- 
els, to fit the blackbody properties of the central star as 
(T,L) = (200,000 K, 2000 Lq). However, uncertainties on the 
photometry, evolutionary models, and extinction law led the 
authors to derive luminosities ranging from 860 - 24,000 L© 
depending on these assumptions. 

Another method for estimating the luminosity of the 
central ionizing source is based the sum of all emitted and 
observed radiation from the nebula. To perform this sum we 
have combined observations from a number of different ob- 
servatories and wavelengths, the results of which are listed 
in Table [21 The lUE observations, which are for a 23 x 10 
arcsec aperture, were scaled up to the whole nebula using 
the intrinsic He ll A1640/ A4686 ratio (-^7 for T ~ 15000 K, 
A^e - 10^ cm"^. Case B, lOsterbrock fc Ferland 200|) and 



the dereddened A4686 line flux from iTsamis et al.l l|2003h . 
The UV - radio continuum emission is primarily free-free 
and bound-free emissi on, the form er of which was estimated 
using the equations ofi AUenl ([1977|), while the latter was cal- 
culated using the 2-photon continuum tool nebcont within 
DlPScQ. The infrared SED was taken from the ISO SWS and 
LWS spectra, with the ISO SWS spectrum was scaled up by 
a factor of 1.45 to match the LWS spectrum. 

The total luminosity of all components is ~5700 L©. 
This is likely to be a lower limit for the luminosity of the 
central star because the nebula and circumstellar disk are 
seen edge on and therefore a fraction of radiation from the 
central star may be escaping along the polar axis. Based on 
the geometry of the model nebula we have used, the open- 
ing angle of the circumstellar disk, and the reprocessing of 
radiation in the bipolar lobes, we estimate that the frac- 
tion of radiation escaping along the polar axis may be as 
much as a factor of two. We therefore estimate the luminos- 
ity of the central star to be ~ 5700 — 11, 400 L©. This value 
is higher than previous estimates, which have either sam- 
pled only one component of the nebula (e.g. lMatsuura et al.l 
I2OO5I . estimated a nebular luminosity of 3300 Lq based on 
infrared dust m odeling) or have bee n estimated from evolu- 
tionary models (jSzvszka et al.ll2009l . estimated 2000 L0). 



3 THE MODEL NEBULA 

The final model described in this paper was built up over 
the course of this work starting with a simple pair of bipo- 
lar lobes and growing in complexity as both the models de- 
manded and more observations were added to the list of 
observations to match. Every part of this model was neces- 
sary in some way to reproduce the observations described 
above. The principle components of the model are a large 
pair of bipolar lobes and a very dense disk-like distribution 
of circumstellar material. In this section we describe the dif- 
ferent components of the NGC 6302 model and how this was 
constructed using the MOCASSIN photoionization code. 

Though observations show the p resence of two pairs o f 
bipolar lobes around NGC 6302 (e.g. iMeaburn et alJIJOOsI '). 
we modeled only one pair for simplicity. Their similar ori- 
entations and the fact that most slits only cover the inner 
regions of the nebula suggest this will not significantly affect 
our results. The presence of large amounts of circumstellar 
material has been inferred from observatio ns of a dark lane 
obscuring the central region of the nebula (jMatsuura et al.l 
|2005|) and ob servations of molec ular material in the core 
of NGC 6302 (|Peretto et al.ll2007l 1. We attempted to model 
this material as both a circumstellar disk or a torus, the 
difference being that a torus has a scale height similar to 
its radial scale length, while a disk is intrinsically fiat with 
a scale length significantly larger than its scale height. We 
found that the observations were better fit by a circumstellar 
disk (this will be discussed further in Section [6. 1|) . 

These structures were all modelled as cylindrically sym- 
metric with the same axes as that of the bipolar lobes. This 
was a necessary simplification for this complicated nebula, 
but it should be noted that radio observations suggest that 
the disk is tilted with respect to the lobes, and m ost-likely 
warped JGomez et al.lll989l : [Matsuura et al.ll2005l '). All these 
components have a number of parameters that were initially 
estimated from previous observations of NGC 6302 but were 
allowed to vary to obtain the best fit during the modeling 
process. The final set of parameters used, and their relative 
uncertainties, are provided in Section [6l 



http://star-www.rl.ac.uk/docs/sun50.htx/sun50.html 



3.1 The bipolar lobes 

The bipolar lobes used in our mode l were developed from 
the model described by iDaval et al.l (|2000| 1 for their model 
of MyCn 18. This geometrical model, shown in Figure [2l is 
based around a truncated paraboloid near the central star 
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Table 2. Total luminosity of NGC 6302 derived from all (dereddened) observed and theoretical components, assuming a distance of 
1.17 kpc llMeaburn et al.ll2008r . See text for a discussion of the individual methods used. 



Component 


Source 


Reference 


Luminosity / Lq 


Ultraviolet line fluxes 


lUE 

ESO 1.52-m 


Feibelman (2001) 


104 


Optical line fluxes 


Tsamis et al. C2003~) 


1324 


Far-red line fluxes 


Siding Springs 2.3-m 


Groves et al. (20021 


73 


Near-IR line fluxes 


UKIRT 


Casassus et al. f2000") 


140 


Free-free continuum 


Theoretical 


Allen figrTI 


122 


Bound-free continuum 


Theoretical 


DIPSO 


837 


Infrared continuum 


ISO 


Molster et al. (2001) 


3090 


Total 






5690 




Figure 3. Visualisation of the three-dimensional density struc- 
ture for the bipolar lobes of NGC 6302 using the model illustrated 
in Figure[2]and using the parameters listed in Table[5] Left: Side- 
on image, showing the hourglass structure. Right: Image of the 
bipolar lobes inclined at an angle of 45° . 



Figure 2. Schematic view of the geometrical structure of an un- 
filled hourglass-shaped bipolar nebula with rotational symmetry 
about the ^-axis. See text for an outline of the parameters govern- 
ing the structure and Table|5]for a list of v alues used in the m odel 
of NGC 6302. Figure adapted from one bv lDaval et~all 1I2OOOI) for 
their model of MyCnlS. 



that develops into a cone-like structure further away from 
the star. The parameters indicated in Figure [2] are the lati- 
tude, 01, of transition from the paraboloid to the cone, the 
latitude, 02, of the end of the cone, the slope, 03, of the 
cone, the semi-height, H , of the hourglass, and the waist- 
radius, Rw, of the paraboloid. The shape of the walls of the 
hourglass can then be described using the equations 



X^ + y^^ f^M. 



R,. 



\tan 



77 55 ] z + Ru 



\z\ < Zp (1) 



where the transition height, Zp, is given by 
. (tan03/tan02) — 1 



z — H- 
^ (tan03/tan0i) — 1 



(3) 



and the radial distance from the hourglass center, R, and 
the latitude, 0, of a point in the nebula are given by 



R = [x + y + z 



2\l/2 



= tan 



(a;2-f y2)i/2 



(4) 



The dimensions of the model nebula were estimated 
from images of NGC 6302 such as that in Fi gure[T]combined 
with an adopted inclination angle of 12.8 ° |Meaburn et al.l 
I2OO5') and distance of 1.17±0.14 kpc (JMeaburnet al .1120081 '). 



We estimated the height of the lobes to be _ff = 7.0 x 10 cm, 
the latitudes of the parabola and cone as 0i = 41.5°, 
02 = 49.6 ° and the slope of the cone as 03 = 66.0 °. These 
quantities were all constrained by images of the nebula and 
were not varied during the modeling process. Estimated 
measurement uncertainties are ±10% for H, ±25% 0i, and 
±10% for 02 and 03. The waist radius, i?^,, was set to be the 
same value as the inner radius of the circumstellar disk and 
could be varied during the modeling process. Using these 
parameters the resulting nebula is quite wide (see Figure [31) 
as is seen in observations of NGC 6302. Finally, the density 
and density gradient in the lobes, as well as whether the 
lobes should be hollow (i.e. only edge-brightened) or filled 
was allowed to vary during the modeling process. 



3.2 The circumstellar disk 

To model the cir cumstellar disk we us ed the flared Keple- 
rian disk model of Pascu cci et al.l (2004). In their model, the 
density distribution of the disk has the form 
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Figure 4. Visualisation of the density distribution in the circum- 
stellar disk seen as a slice through the x-z plane, the image covers 
a region of 1 X 10^^ by 4 X 10^^ cm. The density is represented 
logarithmically in blue and covers the range 10* - 10® cm~^. Note 
the flaring in the outer disk. 



NH(r;z) = No 



exp 



-0.25 TT 2^ 

{zd{r/ri)Sdy 



(5) 



where r is the radial distance and z is the distance from the 
midplane. rd and Zd are disk length and height parameters, 
which determine the scale length of the density distribution, 
and A^o is a characteristic density fo r the disk. The radia l 
density dependence of a = 1 used byJPascucci et al.l (120041 ') 
is also that found bv' Matsuura et al.l ( 20051 ) from a best-fit 
to the spectral energy distribution (SED) in their radiative 
transfer models, though we allowed this to vary in our mod- 
els. The flaring of the disk with distance from the central 
star is given by Zd{r/rd), utilising the flaring parameter 
fd, which controls the opening angle as a function of dis- 
tance from the central star (see Figure |3|. The disk extends 
from an inner radius, Tin, to an outer radius, Tout- AH these 
parameters were allowed to vary in an attempt to find the 
best fitting model. 

Initial parameters f or the circumstel l ar dis k were taken 
from the SED fitting of iMatsuura et ahl (120051 ) such as the 
height parameter, Zd = 5.0 x lO'^^ cm, and the inner radius. 
Tin = 1.0 X 10^® cm . For other parameters we used the initial 
estimates given bv lPascucci et al.l (|2004l ). The outer radius 
was initially set as Tout ~ 1-0 x 10^^ cm. The characteristic 
density, A^o, was also allowed to vary during the modeling 
process. 



3.3 The model grid 

The bipolar lobes and circumstellar disk were mapped onto 
a 2D grid in the x-z plane, utilising the fact that the model 
nebula exhibits both cylindrical symmetry and is symmet- 
ric upon reflection in the z = plane (though, as noted 
above, the observed nebula is not cylindrically symmetric 
and is better described as point-symmetric). The number of 
cells in the grid and the size of the cells was varied during 
the modeling process to ensure that the ionization fronts 
and temperature and density gradients in the model were 
appropriately resolved such that the resulting emission line 
spectrum was a true representation of the model nebula. The 
final model consisted of a 2D grid of 138 x 86 cells with the 
cell width varying from 2 x 10^* in the circumstellar disk, 
to 5 X 10^® in the outer regions of the bipolar lobes. This 
was then mapped onto a 3D grid of 138 x 138 x 86 cells to 



confirm the validity of the 2D models and test the influence 
of density inhomogeneities in 3D (discussed below). 



3.4 The central star 

The central ionizing source of NGC 6302 was mod- 
eled using stellar model atmospheres produced using the 
Tubingen NLTE Mod el Atmosphere Pack aq^ (TMAP, 
iRauch fc DeetienI l2003l : I Werner et all I2003D . While it is 
sometimes true that the assumption of a blackbody is not 
worse than any other assumption in producing a photoion- 
ization model fit, at energies higher than 13.6 eV and 54.4 eV 
the differences between a blackbody flux and a stellar atmo- 
sphere can become so gr eat th at very largo differences will 
resuh (e.g. [Ranch 1997,, l2003l : [Armsdorfer et ah 2002). A 
grid of stellar atmosphere models were tested in our simu- 
lations, with effective temperatures ranging from 50,000 to 
300,000 K and logg' from 5-9. These were tested for both 
solar {[X] = [Y] = [Z] = 0) and metal-poor {[X] = \Y] = 0, 
[Z\ = —1) elemental abundances (w here the solar abun- 
dances used by TMAP are those from lAsplund et al.ll2005l . 
and include all the elements up to nickel) . Model atmosphere 
fluxes were also tried for 'typical' PG 1159 abundance ra- 
tios (He:C:N:0 = 33:50:2:15 by mass), covering the same 
effective temperature and gravity ranges. Given the lack of 
knowledge about the central star of NGC 6302 we allowed 
all these properties to vary in our models. 



3.5 Nebular abundances 

No spatially resolved spectroscopic observations of 
NGC 6302 have been published, so there is no direct 
observational evidence for structures within the nebula with 
different elemental abundances. For this reason we used 
a homogeneous elemental abundance distribution in our 
models, despite mocassin's ability to model regions with 
different chemistries. Our model used 14 elements, including 
all the principle elements that contribute to the heating 
and cooling of the nebula, as well as those responsible 
for the important density and temperature sensitive line 
ratios, and those require d to reproduce the hig h-ionization 
coronal lines observed bv lCasassus et al.l l|2000f) . The initial 
abun dances of He, C, N, O, Ne, S, CI and Ar were taken 
from iTsamis et al.l ((2003), the abu ndances of Mr and A l 
were taken from the analysis of ICasassus et al.l ((20001), 
while a bundances for N a, Si and K were initially set to solar 
values (|Loddersll2003l ). All abundances were initially kept 
constant while structural parameters were being allowed to 
vary, but were then varied to obtain the best fit to the lines 
from each element (see final values in Table [SJ. 



^ Website: http://ast ro.uni-tuebingen.de/rauch/TMAP/TMAP.htmll 
In the framework of the Virtual Observatory (VO, 
http://ww.ivoa.net), all spectral energy distributions from 
the TMAP model grids described here are available in VO 
compliant form from the VO service TheoSSA (http://vo.ari.uni- 
heidelberg.de/ssatr-0.01/TrSpectra.jsp?) provided by the Ger- 
man Astrophysical Virtual Observatory (GAVO, http://www.g- 
vo.org). 
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3.6 Dust modeling 

NG C 6302 is known to have a very large du st component 
(e.g. iMolster et aP I2001I : iKemper et"aDl2002l l that will af- 
fect the radiative transport in the nebula. Any photoion- 
ization model of the nebula must therefore fully consider 
the influence of dust on the ionization structure of the gas. 
The models presented here therefore include a dust com- 
ponent that utilises mocassin's full treatment of dust ra- 
diative transfer (e.g. lErcolano et al.l 12005V The dust model 
will be fully discussed in a future paper, but we outline here 
the essential properties. Dust was included in the form of a 
number of dust species whose s ize distribution was adapted 
from the standard MRN model i|Mathis et al.ll 19771 ). Optical 
constants were obtained from the literature and converted 
to absorption efficiencies using mocassin's light scattering 
codes (that treat both spherical and elliptical dust grains). 
There is limited availability of UV and far-UV optical con- 
stants for common dust species in the liter ature, the mos t 
widely used being the amorphous silicates of lDraind l|2003h . 
which were adopted in our model. Dust absorption within 
the nebula is significant in reducing the ionizing flux, with 
approximately half of all Lyman photons absorbed by the 
dust (in agreement with the observed fraction of luminosity 
coming from the IR SED, see Section [2. 4|) . 



4 THE MODELING PROCESS 

In this section we describe the modeling process that was 
used to constrain the model properties by fitting the dif- 
ferent observations. In particular we describe how different 
model parameters caused variations in the observable lines 
intensities or line ratios and how this led to the final model. 
These dependencies of lines or ratios on certain parame- 
ters were not absolute and changing other parameters could 
sometimes nullify the influence of certain parameters. This 
made the modeling process highly complicated and often 
required the parameter space around a model to be inde- 
pendently investigated at regular intervals in case new pa- 
rameter dependencies had emerged due to changes in the 
model. Furthermore, in a model as complex as this there 
will inevitably be degeneracies in the models that can fit the 
observations, and some lines or ratios may remain not fitted. 
The objective here is to highlight which parts of the model 
are the most constrained by the observations and therefore 
represent the most important features of the model. The fi- 
nal model, while not a perfect fit, avoids many of the poor 
fits that certain parameters resulted in. We have divided this 
discussion, for the purposes of clarity, into the fundamental 
areas of the model (nebular structure, central star proper- 
ties and nebula chemistry), but it should be noted that all 
the parameters were fitted simultaneously. 



4.1 The nebular geometry and density 

The geometry and density of the nebula, specifically its 
two main components, the circumstellar disk and the bipo- 
lar lobes, were refined by fitting the observed density- and 
temperature-sensitive line ratios. These line ratios respond 
to changes in the density and ionization structure of the 



nebula, but also highlight specific regions of the nebula ac- 
cording to their ionization potentials. The diagnostic line 
ratios can therefore be divided into those primarily respon- 
sive to changes in the low-density bipolar lobes, or the high- 
density circumstellar disk. These lists were not independent 
and many line ratios, particularly those sensitive to changes 
in the ionization structure of the nebula, were influenced 
by a large range of parameters. This was further compli- 
cated by the apertures and slits of the different observations 
of NGC 6302, each of which sampled a mixture of the two 
structural components. Despite this, many line ratios very 
highly responsive to specific parameters and by fitting these 
line ratios many structural parameters could be constrained. 
With these restrictions incorporated the variation of the re- 
maining parameters was explored and the best-fitting model 
identified. 

The circumstellar disk is modeled with seven parame- 
ters {No, Tin, Tout, Td, Q, Zd, fd), though the characteristic 
density, A'o, and the inner radius, ri„, were the most infiuen- 
tial in the resulting emission line spectrum and are therefore 
the most well constrained. This is because these parameters 
infiuence the region of the disk exposed to the full ionizing 
fiux, while many of the other parameters affect regions of 
the disk beyond the ionization front. The line ratios that 
were responsive to changes in the circumstellar disk are a 
mixture of those that probe either the ionization structure, 
or regions of high temperature or density. 

The most influential parameter was the characteristic 
density of the circumstellar disk, Nq , which we fit to a value 
of A^o = 80,000 cm~'^, with an accuracy of approximately 
±20%. This parameter was influenced by nearly all the ob- 
served line ratios, though the strongest constraint came from 
the He ll / He I line ratio which was over-estimated by a fac- 
tor 10 if A'^o < 60,000 cm""^, a disparity that could not be 
resolved by adjusting any other parameters. Many of the 
density-sensitive line ratios with high critical densities (i.e. 
those of [CI III], [Ar iv], and [K v] ) could also only be fit 
with a very dense circumstellar disk. Very little variation was 
observed in any of the line ratios when different density pro- 
files were adopted, primarily because the circumstellar disk 
is highly optically thick. The initial choice of Nh on r~^ for 
the circumstellar disk radi al density profile was m ade on the 
basis of the SED fitting of iMatsuura et al.l l|2005f ). but r"^/^ 
or r~^ density profiles led to little difference. 

The inner radius of the circumstellar disk, ri„, was also 
very influential as it dictates the strength of the ionizing 
flux at the circumstellar disk. It was best constrained by a 
mixture of the [O ll] (3726 + 3729) / (7320 + 7330) and 
[Ar iv] density-sensitive line ratios. The former was over- 
estimated by a factor 5 if ri„ > 1.4 x 10^^ cm, and the latter 
increased significantly as rin was decreased. The final value 
of Tin = 1.2 X 10^^ cm was chosen to balance these two ratios 
though it could be decreased if the [Ar iv] line ratio could 
be reduced by other methods. At a distance of 1.17 kpc this 
is equivalent to an angular radius of 0.68". 

The thickness of the circumstellar disk (determined by 
Zd) infiuenced a large number of line ratios, most critically 
the [O ll] 3729 / 3726 ratio which became irretrievably over- 
estimated if the disk thickness was increased, and the [N ll] 
5755 / (6584 -I- 6548) ratio that greatly increased if the 
thickness was significantly decreased. The best fit value of 
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Zd ~ 6.0 X 10^^ cm was constrained to an uncertainty of 
±50%. 

The other disk parameters only influence material be- 
yond the ionization front and therefore had a lot less in- 
fluence on the modeled spectrum. The flaring parameter, 
fd, was only constrained to be J5 f .f by the [O in] / [O ll] 
and [S ll] 4068 / (6731 -I- 6717) line ratios (which both in- 
creased significantly if the parameter was reduced), but no 
upper constraints were identified. It was therefore kept at 
its initial value, fd = 1.125. The length parameter, r^, and 
the outer radius, rout, influenced a small number of line ra- 
tios, but none sufficiently enough to constrain them. The 
former parameter was therefore kept at its initial value of 
Td ~ 5.0 X lO^'' cm, while the latter parameter was adjusted 
based only on the results of radiative transfer dust-fitting 
(not discussed in this paper) and therefore we list its value 
as unconstrained here. 

The circumstellar material was modeled here as a flat 
disk, but attempts were made to model it as a torus 
(i.e. a circumstellar distributio n with a greater height and 
smaller length), as s uggested bv lMatsuura et al.l ((20051) and 
IPeretto et al.l (|2007|) . However, increasing the disk thickness 
to that suggested bv lPeretto et al.l ([20071) resulted in an un- 
flttable [O ll] 3729 / 3726 ratio (as noted above). We there- 
fore suggest that certainly in the inner regions where the 
ionized gas traces the structure of the circumstellar mate- 
rial, it is in the form of a thin disk. In the outer, neutral 
regions the circumstellar material may better represent a 
torus if traced by neutral material or dust. 

The shape and dimensions of the bipolar lobes were 
adopted from images of the nebula s uch as that in Figure [l] 
with the distance to the nebula from lMeaburn et al.l (|200g ). 
These parameters were flxed throughout the modeling pro- 
cess and only the density of the lobes and the density pro- 
file were varied. The lobes were initially modeled with a low 
density of 2000 cm~^ constrained by a number of line ratios 
that are particularly sensitive to density such as the [O ll] 
3729 / 3726, [O iii] 4363 / 5007, [O iii] 52 / 88 /^m, and [S ii] 
(A4068 -f A4076) / (A6731 4- A6717) ratios (aU sensitive to 
densities < 5000 cm""^). However, none of these line ratios 
offered a particularly strong constraint on the final density, 
with the [O ll] 3729 / 3726 line ratio requiring densities in 
excess of 1000 cm"^, while the [O ill] 52 / 88 /^m line ra- 
tio increased if the density was significantly increased above 
4000 cm~^. A final uncertainty of ±50% is estimated for the 
final value of 2000 cm~^ based on these constraints. No con- 
straint could be found on the density proflle of the bipolar 
lobes so they were kept at a constant density. 

In response to an initially poor fit to a number of 
density-sensitive line ratios such as those of [Ar iv] and 
[K v], as well as the temperature-sensitive [N ll] line ratio, 
a third component was added to our model. This region, 
dubbed 'the outfiow', is an inner region of the bipolar lobes 
with a higher density, N out flow, and with the same geometry 
but with an outer radius, rout flow, at which the density of 
the bipolar lobes dropped to the standard value. Introduc- 
ing this component greatly reduced the [N ll], [Ar iv], and 
[K v] , line ratios that had previously been significantly over- 
estimated. The final values of the two outfiow parameters, 
Noutfiou, = 20, 000 cm-^ and rout flow = 2.0 x lO" cm (11.5" 
at 1.17 kpc), are well constrained to ±20%, due mainly to 
the first two line ratios, while the third offered less of a con- 



straint. As an example of the modeling process Table |3] lists 
the line ratios most responsive to the introduction of, and 
changes to, the outflow component. 

4.2 The central star and the IR coronal lines 

Without any clear observations of the central star its pa- 
rameters were particularly unconstrained at the beginning 
of the modeling process. While there had been many sug- 
gestions in the literature f or a very high central star temper- 
ature (e.g. I Ashley fc Hvland 1988 : Casassus et al. 20001 ) we 
began the modeling process with a typical FN central star 
with a blackbody temperature of 90,000 K and solar abun- 
dances. While a large number of line ratios were responsive 
to changes in the central star temperature, the ionization 
structure diagnostics He ll A4686 / He I A5876 and [O ill] 
(A4959 ± A5007) / (A3726 ± A3729) were particularly re- 
sponsive to it and would not come close to being flt unless 
the blackbody effective temperature was increased to at least 
~150,000 K. 

At this temperature, with the peak of the SED in the 
far-UV, further changes produced little response in any line 
ratios. With the exception of the total H/3 luminosity of the 
nebula, the stellar luminosity was also very unconstrained at 
these temperatures. The final stellar luminosity was there- 
fore based on the best fit to the H/3 line strength. For many 
FNe this may be a relatively unsatisfactory result, but given 
the high precision to which the distance to NGC 6302 is 
known, its absolute H/3 luminosity and therefore the lumi- 
nosity of its central source, is quite well constrained. 

To constrain the stellar temperature further we used ob- 
servations of a group of emission lines from high-ionization 
stage ions (the 'infrared coronal lines') t hat have been ob- 
served in the centre of NGC 6302 (e.g. Ashlev fc Hvlandl 



1 19881 : iFottasch et"aLlll996l : ICasassus et al.ll200Gl ). This group 
includes emission from ionization stages ranging up to Si®^ 
and are unique to a very small number of PNe, with 
NGC 6302 showing the highest ionization species ever found 
in a FN. These were originally susp ected to originate in high- 
velocity shocks in the n ebula (e.g. iMeaburn fc Walshlll98Gl : 
iLame fc Ferlandl ll99H ), but p hysical argument s against 
such levels of shock excitation (jOliva et al.lll996l ) and the 
lack of significantly broad wings on many of the lines 
(jCasassus et al.ll200(l ) argues for a photoionized origin. 

The ionization potentials of many of these species lie 
in the extreme UV or X-ray region (see Table |4)| and there- 
fore place strong restrictions on the form of the high-energy 
ionizing flux from the central star. As discussed above, the 
ionizing flux from the central star was modeled using NLTE 
stellar model atmosphere fluxes for a range a range of cen- 
tral star abundances and surface temperatures. The effects 
of different temperatures and abundances on the fluxes of 
the six highest ionization stage ions is shown in Table U 
The effects of these changes on the lower ionization-stage 
nebula diagnostics was negligible. 

Models using either a solar abundance or metal-poor 
stellar atmosphere were able to produce fits to many of the 
lines but were unable to reproduce the lines from the highest 
ionization stages, particularly the [Si ix] 3.93 /mi line. The 
strength of this line is underestimated by a factor ^ 10'' even 
at Teff — 250, 000 K and an order-of-magnitude fit could 
not be obtained for any models with Te// < 400, 000 K. At 
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Table 3. The effects of introducing and varying the parameters of the outflow component, on the most responsive density and temperature 
sensitive line ratios. Note how the majority of line ratios decrease as the outflow is introduced and its density increased, but then increase 
at very high outflow densities. 



Line ratio 



Observed 



Models with outflow density: (cm '^) 



[N II] 5755 / (6585+6546) 

[O III] (4959+5007) / [O ii] (3726+3729) 

[S II] 4068 / (6731 + 6717) 

[Ar iv] 4741 / 4711 

[K v] 4163 / 4123 





No outflow 


10,000 


20,000 


40,000 


0.0290 


0.133 


0.0227 


0.0231 


0.0278 


39.1 


401 


259 


87.8 


60.6 


0.485 


2.16 


1.13 


0.777 


1.38 


2.01 


7.32 


4.96 


4.29 


5.10 


1.07 


8.38 


4.61 


4.32 


4.12 



Table 4. Fits to the observed infrared coronal lines from NGC 6302 using NLTE stellar model atmospheres with solar or hydrogen- 
deficient abundance compositions (see Section 13.411 . Results for three effective temperatures are shown for each stellar abundance model 
(both with logq = 7). Line fluxe s listed as 0.0 are at least a factor of 10'* smaller than the observed flux. Dereddened line fluxes are 
taken from lCasassus et al.l 1120001 ). All line fluxes are given in units of 10^^^ erg cm~^ s~^. Ionization potentials (IP) are also listed for 
each prior ionization stage. 



Line 


IP 

(eV) 


Observed 


Models 




Solar composition 'Typical' H-def composition 




200kK 220kK 240kK 200kK 220kK 240kK 



[Mg v] 5.60 /im 
[Mg vil] 5.51 lira 
[Mg viii] 3.03 Aim 
[Si vi] 1.96 /im 
[Si Vll] 2.47 /im 
[Si ix] 3.93 Atm 



109 
187 
225 
167 
205 
303 



16.9 
13.0 
1.81 
23.5 
20.5 
0.030 



6.33 
0.67 
0.0 
14.1 
0.0 
0.0 



21.0 
1.52 
0.0 
21.7 
0.13 
0.0 



222 
28.1 
0.32 

429 
4.51 
0.0 



7.4 
9.42 
0.91 
31.6 
4.23 
0.002 



26.6 
9.58 
2.82 
73.4 
32.1 
1.04 



0.23 
0.58 
4.67 
1.34 
9.32 
139 



such high stellar temperatures the lower ionization coronal 
lines (e.g. [Mg v] 5.60 /im and [Si vi] 1.96 nm) become over- 
estimated by factors of several 100 or more. Models using 
hydrogen-deficient stellar atmospheres (Table |4)| produced 
substantially better fits to these lines, with the best fit ob- 
tained using a 220,000 K stellar atmosphere. The fits to the 
magnesium lines presented in Table |4] are very good. The 
fits to the nebular silicon lines are slightly worse, with the 
2.47 Atm line fitted well, but the 1.96 /im line is overestimated 
by a factor of three and the 3.93 ^^m line is overestimated 
by a factor of thirty. 

FigureOshows the high-energy SEDs of 220,000 K solar- 
abundance and H-deficient model atmospheres. At the en- 
ergies of the Mg^^ and Si"*"^ ionization potentials, there is 
a significantly higher flux in the H-deficient model atmo- 
spheres compared to the solar abundance model due to the 
lack of opacity from hydrogen. The absence of a number of 
low-abundance metals in these model atmospheres results in 
a reduced atmospheric opacity and therefore a higher flux 
at energies > 250 eV. Including them will likely reduce the 
photon flux at 330 eV and therefore resolve the problem of 
the over-predicted flux in the [Si ix] 3.93 /.im line in our 
models. 

With the exception of a blackbody (which provides a 
better flt to the IR coronal lines, but i s not an accurate 
repre sentation of the ionizing flux, see e.g. lArmsdorfer et al.l 
[2003), the 220,000 K H-deficient stellar atmosphere provides 
the best fit to the observed IR coronal lines. Models with dif- 
ferent surface gravities were tested but produced very minor 
changes in the model line fluxes. A value of log g — 7.0 pro- 
duced the best flt, but is poorly constrained. The fit to the 



firP^r'T' 




Solar abundance model atmosphere. log(g} = 7, T=220.000 K 
H-deficient model atmosphere, log(g) = 7, 1=230,000 K 
Blackbody, T=220,000 K 



50 



100 



200 



hv [eV] 

Figure 5. 220,000 K solar abundance and hydrogen-deficient 
model atmosphere SEDs compared with a blackbody at the same 
temperature. The frequencies corresponding to the ionization po- 
tentials of Mg+^ and Si"*"^ are shown. 



higher ionization stages of the silicon ions is not ideal, but it 
provides a closer fit than any of the solar composition mod- 
els and is strong evidence that the central star of NGC 6302 
is H-deficient. 
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4.3 The nebular chemistry 

Once the main nebular and central star parameters were de- 
termined by fitting the density- and temperature-sensitive 
line ratios, attempts were made to achieve good fits to 
the strengths of individual emission lines by altering the 
chemical abundances of the nebula, which had originally 
been set to li t eratu re va lues for NGC 6 3 02, p rimarily from 
iTsamis et all (|2003l ) and lCasassus et al.1 l|2000l ). The major- 
ity of abundances required small changes to fit the emission 
line fluxes, resulting in moderate variations in the nebula 
temperature structure which were then resolved by making 
other changes to the nebula structure. The sodium, mag- 
nesium, aluminium and silicon abundances were unchanged 
from their initial values either due to satisfactory fits, or 
insufficient numbers of lines to fully constrain their abun- 
dances. The final nebula abundances are listed in Table (5] 

A full discussion of the dust modeling process will be 
reserved for a future paper, but as noted earlier dust was 
included in these models. This caused variations in a large 
number of line ratios due to significant changes in the neb- 
ular ionization structure. 



5 MODEL RESULTS 

In this section we present the best fitting model of NGC 6302 
and compare the predicted emission-line fiuxes from MO- 
CASSIN with those observed. The list of parameters used in 
the final model outlined in Section 13.11 of NGC 6302 are 
provided in Table [S] We also note how well constrained we 
estimate each parameter to be. The uncertainty on each pa- 
rameter is estimated as the degree to which the parameter 
could vary without adversely affecting the model. For neb- 
ular abundances the uncertainty is based on the differences 
between observed and predicted line strengths for each ele- 
ment. 



(UKIRT). We note a significantly higher discrepancy be- 
tween model and observations for line fluxes measured from 
ISO spectra than for line fluxes measured from the UKIRT 
spectra. This could be due to the combined effects of a 
misaligned pointing used for the ISO observations and the 
relatively small size of the ISO S WS aperture, as in th e 
case of the NGC 391 8 mod eled bv lErcolano et al.l (|2003h . 
iBeintema fc PottaschI (|l999i ) noted a significant disconti- 
nuity between short and long wavelength sections of the 
SWS exposures, which they attributed to misalignment of 
the pointing. Combined with the 3" pointing error on the 
14 X 20 arcsec ISO SWS aperture positioned over the central 
10" diameter peak emission region suggests errors arising 
from such an offset could be high. As with all our predicted 
line fluxes, we have reproduced the relevant aperture or slit 
used for the observations and therefore any inaccuracies in 
the observational pointing would seriously affect a compar- 
ison of line fluxes, especially for an object like NGC 6302 
with such strong density and temperature contrasts in its 
structure. The 750 LWS aperture is much larger, making 
any potential pointing error smaller, and we find a smaller 
discrepancy between our model results than for the SWS 
measurements, supporting this. 

Many infrared lines from neutral species are predomi- 
nantly emitted from the photo-dissociation regions (PDRs) 
of nebulae fe.g. lLiu et al.ll200ll : IVasta et al.ll2010 ) which are 
not treated by the current version of MOCASSIN. Despite this 
Table |6] indicates that the [O l] 63 /im and 145 /im as well 
as the [C ll] 158 pim lines show good agreement between the 
predicted and observed values, despite the fact that they 
originate from the neutral FDR. The current version of MO- 
CASSIN therefore appears to be matching these lines ade- 
quately. To complement this and for comparison with fu- 
ture studies we use our model to predict fluxes for two FDR 
lines without published fluxes, but which may be measured 
by future instruments. For the [C l] 369 and 609 /im lines 
we derive fluxes for the entire nebula of 7.80 x 10^^^ and 
1.45 X 10"^^ ergs cm"^ s"^ 



5.1 Fits to the emission line spectrum 

Predicted emission line fiuxes for our best-fitting model are 
given in Table [H] with strengths given relative to the in- 
trinsic dereddened H/? flux where possible. The majority 
of the line strengths presented are in reasonable agreement 
with the observations, with fits to within 30%. In pa rticular, 
nearly all the UV and optical line fluxes from Tsa mis et al.l 
1(2003) ^"^^ "^^11 reproduced. We obtain a reasonable fit to 
the N V 1240 A resonance line which is often overestimated 
in photoionization modeling. Resonance lines must travel a 
much greater distance to escape from the nebula than do 
forbidden lines because of the large number of scatterings 
that they undergo and therefore can suffer greater atten- 
uation due to dust. MOCASSIN treats these lines separately 
and their absorption by dust is accounted for in this model 
l|Ercolano et al.fl2005l ). 

The majority of infrared fine-structure line flux mea- 
surements are taken fro m spectra taken with the ISO SWS 
and LWS , prese nted by iBeintema fc FottaschI (|l999l ) and 
ILIu et al.l (|200ll ). respectively. The remaining infrared lin e 
flux measurements were presented bv lCasassus et al. ' (2000), 
from observations made with the UK Infrared Telescope 



5.2 Fits to the density- and temperature-sensitive 
line ratios 

Table [7] lists the model emission line ratios that are sensi- 
tive to changes in the density, temperature, or ionization 
structure of the nebula, and compares them to the observed 
line ratios of NGC 6302. Despite the wide range of elec- 
tron densities and temperatures noted in the literature for 
NGC 6302 we have been able to reproduce the majority of 
emission line ratios using our model. The flts generally show 
good agreement to within 20%, with a small number only 
fit to within a factor of two. The most notably bad fit is the 
density-sensitive [K v] 4163 / 4122 line ratio, which is over- 
predicted by a factor of ^^4 and proved particularly difficult 
to fit. 

The [Ar iv] and [K v] line ratios are both over-predicted 
by our model suggesting we have over-estimated the den- 
sity of the circumstellar disk. However, many other line ra- 
tios constrained the density of the circumstellar disk and 
no suitable combination of parameters could be found us- 
ing a lower-density circumstellar disk. Increasing the inner 
radius of the circumstellar disk offered a potential solution 
to fit these two line ratios, but increasing this parameter to 
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Table 5. Model parameters for the best fitting final model of NGC 6302 with the origin of the value and the uncertainty on the value 
(the listed uncertainties do not take into account the uncertainty of ~10% on the distance to NGC 6302). All nebula abundances were 
constrained by the model fits, but we also note the values that remained unchanged from the initial values listed in Section l3.5l Elemental 
abundances are listed as [X/H] , with enhancements or depletions noted in addition to the original source. 



Model parameter 


Value 


Origin 




Uncertainty 


Bipolar lobes: 










Height of lobes, H 


6.8 X 10" cm 


Estimated from 


images 


±10% 


Transition latitude, cp\ 


41.5° 


Estimated from 


images 


±25% 


End of cone latitude, (/)2 


49.6° 


Estimated from 


images 


±10% 


Slope of cone, 03 


66.0° 


Estimated from 


images 


±10% 


Transition height, Zp 


4.0 X 10^^ cm 


Estimated from 


images 


±25% 


Waist radius, R^ 


1.2 X 10^6 cm 


Model fit 




±20% 


Lobe density 


2000 cm-3 


Model fit 




±50% 


Density profile 


constant 


Model fit 




- 


Outflow density 


20,000 cm-3 


Model fit 




±20% 


Outflow outer radius, Voutflow 


2.0 X 10^^ cm 


Model fit 




±20% 


CircumstoUar disk: 










Inner radius. Tin 


1.2 X lOl'' cm 


Model fit 




±20% 


Outer radius. Taut 


3.0 X lOl'^ cm 


Not well constrained 


- 


Length parameter, r^ 


5.0 X lO^'^ cm 


Model fit 




±50% 


Height parameter, z^ 


6.0 X 10^5 cm 


Model fit 




±20% 


Flaring parameter, /^ 


1.125 


Model fit 




±50% 


Characteristic density, A^o 


80,000 cm-3 


Model fit 




±10% 


Nebula abundances: 










Hydrogen 


1.0 


- 




- 


Helium 


0.153 


Model fit 




±10% 


Carbon 


2.2 X 10-* 


Model fit 




±25% 


Nitrogen 


3.9 X 10-* 


Model fit 




±25% 


Oxygen 


5.1 X 10"* 


Model fit 




±25% 


Neon 


2.4 X 10~* 


Model fit 




±25% 


Sodium 


2.1 X 10-*^ 


Model fit 




±25% 


Magnesium 


1.9 X 10-^ 


Model fit 




±25% 


Aluminium 


3.0 X 10~* 


Model fit 




±50% 


Silicon 


3.6 X 10-^ 


Model fit 




±50% 


Sulphur 


2.5 X 10-^ 


Model fit 




±25% 


Chlorine 


2.5 X 10-^ 


Model fit 




±10% 


Argon 


1.2 X 10~^ 


Model fit 




±25% 


Potassium 


5.0 X 10-^ 


Model fit 




±25% 


Central star properties: 










Abundances 


He:C:N:0 = 33:50:2:15 


Model fit 




- 


Effective temperature 


220,000 K 


Model fit 




±10% 


Luminosity 


14,300 Lq 


Model fit 




±20% 


logs 


7.0 


Model fit 




±20% 



Tin ~ 1.5 X 10^^ cm resulted in extreme values of the [O ll] 
(3726 + 3729) / (7320 + 7330) line ratio. The introduction 
of an 'outflow' component (see Section |4Tj went a long way 
to resolving these discrepancies, but no simple model could 
be found to fit these lines accurately. Given this situation 
we were forced to leave these line ratios not fitted by our 
final model, but Section 15.2.11 discusses a possible solution 
to this problem. The poor fit to the He ll / He I line ratio 
is also worrisome, but again Section [5. 2. II presents a way to 
resolve this problem. 



The majority of other line ratios are fit well and in some 
cases offer strong constraints on the geometry or density 
distribution of the nebula (see Section |4]T] for a discussion). 



5.2.1 The effects of density inhomogeneities on the 
emission line spectrum 

A number of the density- and temperature-sensitive line ra- 
tios discussed in Section 15.21 such as the He ii 4686 / He i 
5876, [Ar iv 4741 / 4711, and [K v] 4163 / 4122 line ra- 
tios, are not well fit by our model. Exploring the parameter 
space in our model (see discussion in Section [4. If) revealed 
that almost every line ratio could be fit with some realis- 
tic combination of model parameters, but that no simple 
solution could be found that would fit all the line ratios. 
This does not necessarily suggest that a purely photoioniza- 
tion model solution does not exist for NGC 6302. A possible 
solution to this problem is that the morphology and den- 
sity distribution of our model is an over-simplification of 
the true density distribution. The inhomogeneous appear- 
ance of the bipolar lobes of NGC 6302 (e.g. Figure [l| sug- 



© 0000 RAS, MNRAS 000, 000-000 



12 N.J. Wright et al. 

Table 6. Observed and predicted emission line fluxes for NGC 6302. The H/9 line flux has been dereddened and modelled using a 
distance of 1.17 kpc. The UV and optical emission lines are listed in units where I(H/3) = 1 while the IR lines longward of 0.9 /xni are 
given in units of 10~^^ erg cm~'^ s~^. See Section 12.31 for a summary of the observations being matched here and the source of each 
observation. All relevant slits and ape rtures were simulated in MOCA SSIN to accurately repr o duce the neb ular areas from which the lines 
were obse r ved. R eferenc es are: R94 =lRowlands et al.l lll994l'). B99 =lBeintema k. Pottasciil l|l999l 'l. COO = ICasassus et al.l l|200d '). LOl = 
iLiu et al.1 | |200J) , G02 = lGroves et al.l ||2002| '). and T03 = lTsamis et al.l ||2003| '). 



Line 



Observed Model Reference Line 



Observed 



Model 



Reference 



erg cm 



H^/10 
H/3 4861 
He I 5876 
He I 4471 
He 11 4686 
[C ii] 157.7 ^lm 
[N i] 5199 
[N ii] 5755 
[N ii] 6548 
[N II] 6584 
[N III] 57.3 ^lra 
N iv] 1486 
N V 1240 
[O i] 5577 

6300 

63.2 ^m 

145.5 /im 

3726 

3729 

7320 

7330 
51.8 /im 
S.4 /im 
5007 
4959 
4363 



[Oi] 
[Oi] 
[Oi] 
[Oil] 
[On] 
[On] 
[On] 
[O III] 
[O III] 
[O III] 
[O III] 
[O in] 



[O iv] 25.9 Aim 



[Nc ii] 
[Nc in] 
[Nc III] 
[Ne III] 
[Ne IV] 
[Ne IV] 
[Ne IV] 



12.8 ^m 
15.6 /im 
36.0 /^m 
3967 
2423 
4724 
4725 



776 


852 


T03 


1.00 


1.00 


- 


0.175 


0.144 


T03 


0.0592 


0.0489 


T03 


0.753 


1.18 


T03 


13.7 


22.7 


LOl 


0.108 


0.102 


T03 


0.201 


0.205 


T03 


1.74 


2.19 


T03 


5.22 


6.67 


T03 


167 


248 


LOl 


5.34 


5.86 


T03 


13.82 


19.0 


T03 


0.00232 


0.00280 


T03 


0.243 


0.338 


T03 


284 


135 


LOl 


10.5 


8.27 


LOl 


0.300 


0.158 


T03 


0.140 


0.0517 


T03 


0.0857 


0.0866 


T03 


0.0728 


0.0713 


T03 


157 


177 


LOl 


49 


20.2 


LOl 


12.9 


13.8 


TOl 


4.29 


4.62 


T03 


0.396 


0.316 


T03 


322 


472 


B99 


25.6 


136 


COO 


377 


186 


B99 


21.0 


60.1 


B99 


0.276 


0.338 


T03 


7.95 


4.96 


T03 


0.0208 


0.0177 


T03 


0.0182 


0.0150 


T03 



[Ne v] 3426 
[Ne v] 14.3 fira 
[Ne v] 24.3 ^tm 
[Ne vi] 7.64 fj.m 
[Na vi] 8.64 /^m 
[Na vii] 4.69 /^m 
[Mg v] 5.60 ^lra 
[Mg Vll] 5.51 ^na 
[Mg viii] 3.03 Aim 
[Al v] 2.88 /im 
[Al vi] 3.66 /im 
[Al viii] 3.69 ^lm 
[Si ii] 34.6 ^lra 
[Si vi] 1.96 ^lm 
[Si vii] 2.47 ^lra 
[Si ix] 3.93 ^lm 
[S ii] 4068 
[S II] 6716 
[S II] 6731 
[S III] 6312 
[S III] 9533 
[S III] 18.7 Atm 
[S IV] 10.5 ^lm 
[CI III] 5517 
[CI III] 5537 
[Ar ii] 6.98 Atm 
[Ar III] 7751 
[Ar III] 7135 
[Ar iv] 4741 
[Ar iv] 4711 
[Ar Vl] 4.53 fira 
[K v] 4123 
[K v] 4163 
[K vii] 3.19 Aim 



2.33 


6.12 


T03 


634 


1585 


B99 


308 


434 


B99 


315 


1478 


COO 


10.7 


7.98 


B99 


4.00 


2.13 


COO 


16.9 


26.6 


COO 


13.0 


9.58 


COO 


1.81 


2.82 


COO 


0.0124 


0.0742 


COO 


0.141 


0.102 


COO 


0.16 


0.0649 


B99 


17.4 


325 


B99 


23.5 


73.4 


COO 


20.5 


32.1 


COO 


0.030 


1.74 


COO 


0.145 


0.218 


T03 


0.101 


0.0893 


T03 


0.198 


0.191 


T03 


0.0584 


0.0328 


T03 


0.717 


0.666 


B99 


57.7 


55.3 


B99 


44.0 


18.9 


COO 


0.00309 


0.00391 


T03 


0.00861 


0.0127 


T03 


51.0 


45.8 


B99 


0.0738 


0.0474 


G02 


0.237 


0.198 


T03 


0.209 


0.191 


T03 


0.126 


0.0445 


T03 


29.0 


60.3 


COO 


0.00382 


0.00150 


T03 


0.00407 


0.00648 


T03 


0.279 


0.323 


COO 



gests that the density distribution is unlikely to be smooth, 
with many clumps, knots and small-scale structures visible 
fe.g. iMatsuura et al.ll2005l . and Figure [1}. Density enhance- 
ments can cause variations in density-sensitive line ratios 
as well as temperature-sensitive ratios though self-shielding. 
High densities can also induce coUisional suppression of some 
coUisionally-excited line, potentially invalidating the tem- 
perature sensitivity of some line ratios. The different densi- 
ties determined from the [O ill] far-IR line ratio compared 
to the [CI hi] and [Ar iv] line ratios, despite their similar 
ionization potentials (35.1eV compared to 23.8 and 40.7eV 
respectively), also suggests that moderate density inhomo- 
geneities exist in the nebula and are responsible for quench- 
ing of the [O hi] far-IR l ines by collision al de-excitation in 
the high-density regions (jLiu et al.ll200H ). 

Attempting to model this structure in detail would be 
impractical and largely unproductive compared to our at- 
tempts to characterise the overall structure of the nebula in 
Section 14.11 However, exploring the influence of such inho- 
mogeneities will allow us to test the theory that they are 



responsible for the discrepancies between our model and the 
observations. To simulate these density inhomogeneities a 
number of 'knots' of various sizes and densities were simu- 
lated in our model using mocassin's subgrid feature to re- 
solve small-scale structures. For simplicity these knots were 
simulated as spherical over-densities in the structure of ei- 
ther the bipolar lobes or the circumstellar disk. Two parame- 
ters were used to model these structures, fp, a density factor 
representing the peak over-density of the knot , and r^not , the 
radius of the knot. All the knots were modeled using a r~^ 
radial density profile. These parameters were allowed to vary 
over the ranges 1.0 < fp < 3.0 and 10^® < rknot/cm < lO^'^, 
the latter of which was estimated from images of the neb- 
ula. The knots were distributed randomly across either the 
bipolar lobes, the circumstellar disk, or both, and the total 
number of knots was varied from 5-50. 

The effect of these density enhancements on the neb- 
ula ionization structure depends on their position in the 
nebula and their density, but their main effect is to cre- 
ate shielded regions where lower ionization stages can ex- 
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Table 7. Comparison of model and observed ratios (with uncertainties) for density-, temperature-, and io nization-sensitive line ratios 



fo r NGC 6302. Also l isted are the quoted densities a nd tem peratures from the literatu re. References : D73 = iDanziger et al. | |l97i ), B8 2 
= lBarraI etHl l|l982h . R94 = [Rowlands et al.l l|l994l '). LOl = lLoidl et all HfloJ), G02 = lGroves et~all l|2002l ). T03 = lTsamis et al.l |20o3 ). 



Line ratio 



Observed 



Model 



Density (cm •^) 



Reference 



[C III] 1907 / 1909 

[O II] 3729 / 3726 

[O II] (3726 + 3729) / (7320 

[O III] 52 / 88 ^m 

[Ne v] 24.3 / 14.3 /im 

[S II] 6731 / 6717 

[CI III] 5537 / 5517 

[Ar iv] 4741 / 4711 

[K v] 4163 / 4122 



7330) 



1.1 ±0.2 
0.449 ± 0.05 
2.78 ±0.08 
3.20 ±0.12 
0.49 ±0.05 
1.96 ±0.14 
2.54 ±0.28 
2.01 ±0.18 
1.07 ±0.14 



0.522 
0.327 
1.33 
8.76 
0.27 
2.14 
3.25 
4.29 
4.32 



16000 
5000 
5750 
1380 
10000 
12900 
22450 
14900 
10000 



B82 
G02 
T03 
LOl 
B99 
T03 
T03 
T03 
T03 



Line ratio 



Observed 



Model Temperature (K) Reference 



[N ii] 5755 / (6584 + 6548) 

[O i] 5577 / 6300 

[O III] (4959 + 5007) / 4363 

[O III] 4363 / 5007 

[S II] 4068 / (6731 + 6716) 



0.0290 ± 0.005 

0.0126 ± 0.0020 

45.7 ±6.5 

0.029 ±0.010 

0.485 ± 0.082 



0.0231 
0.0083 

58.3 
0.0229 

0.777 



14225 
10000 
18400 
17400 
10000 



T03 
G02 
T03 
D73 
T03 



[S III] 6312 / 9533 








0.049 ±0.14 


0.0475 


18300 


B99 


Line ratio 








Observed 


Model 




Reference 


He I 5876 / 4471 

He II 4686 / He I 5876 

[O III] (4959 + 5007) / [O 


II] (3726 


+ 


3729) 


2.96 ±0.21 
4.30 ±0.30 
39.1 ±6.4 


2.94 
8.19 
87.8 




T03 
T03 
T03 



ist. The most notable effects are therefore on the line ratios 
that probe the ionization structure of the nebula, particu- 
larly the He ll / He I line ratio which is significantly re- 
duced if knots are introduced into the circumstellar disk. 
The other ionization-dependent line ratio, [O ill] (4959 + 
5007) / [O ll] (3726 + 3729), was also found to decrease if 
knots were added to either the lobes or the disk, which had 
the effect of worsening the fit to this line ratio. However the 
magnitude of this change was smaller than that of the he- 
lium line ratio, suggesting it could be countered with other 
model changes. 

Changes to density-sensitive line ratios were also ob- 
served when knots were added to the bipolar lobes, particu- 
larly the [CI III] , [Ar iv] , and [K v] fine ratios that aU probe 
high density regions, which were found to move closer to 
the observed ratios. The [O ll] and [S ll] line ratios, sensitive 
to lower-density regions, were less responsive. Adding knots 
to the circumstellar disk did not induce significant changes 
in these density-sensitive line ratios, most probably because 
either the upper critical densities of these line ratios are 
lower than the typical densities in the circumstellar disk or 
because the majority of the circumstellar disk is in a low 
ionization state. 



Table 8. Integrated ionic abundance ratios for He, C, N, O, Ne, 
S, CI, and Ar, derived from model i onic fractions and compared 
to the results of lTsamis et al.l ||2003| '). 



Ionic ratio 



He+ / H+ 
He2+ / H+ 
C2+ / H+ 
C3+ / H+ 
N+ /H+ 
N2+ / H+ 
N3+ / H+ 
N4+ / H+ 
0+ /H+ 
02+ / H+ 
0=*+ / H+ 
Ne2+ / H+ 
Ne3+ / H+ 
S+ /H+ 
S2+ / H+ 
Cl2+ / H+ 
Ar2+ / H+ 
Ar3+ / H+ 



Model 



Observed 



0.0645 


0.0658 


0.0477 


0.0696 


1.82(-5) 


3.08(-5) 


9.51(-5) 


1.43(-5) 


1.34(-4) 


4.48(-5) 


7.44(-5) 


1.03 (-4) 


3.75(-5) 


7.80(-5) 


4.06(-5) 


9.45(-5) 


3.00(-4) 


1.67(-4) 


4.17(-5) 


9.35(-5) 


1.88(-5) 


8.22(-5) 


1.19(-5) 


2.65(-5) 


5.04(-5) 


3.53(-5) 


1.28(-6) 


8.20(-7) 


1.98(-6) 


1.57(-6) 


1.13(-7) 


2.76(-8) 


7.10(-7) 


6.51(-7) 


3.39(-7) 


7.75(-7) 



In summary, introducing small-scale density inhomo- 
geneities or 'knots' in the nebula structure could resolve 
some of the remaining discrepancies between our model and 
the observations. Knots in the circumstellar disk could re- 
solve the over-predicted ionization-sensitive helium line ra- 
tio, while knots in the bipolar lobes could resolve differences 
in line ratios sensitive to high-densities. Despite this, a num- 
ber of the remaining line ratios could not be matched within 



the parameter space explored (e.g. the low-density-sensitive 
[O II] (3726 + 3729) / (7320 + 7330) line ratio). 



5.3 The nebular ionization structure 

A comparison of the relative ionic abund ances in our model 
with those determined by iTsamis et al.l l|2003r ) is shown in 
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Table [HI The ionization structure derived by those authors 
is in good agreement with our model, though there are a 
number of significant diff'erences. The ionic fractions of O "*" 
differ by a factor of four, while the fraction of Cl^^ in our 
models is a factor of five larger than the empirical observa- 
tional value. Fits to the [O iv] and [CI lll] lines in Table [6] 
are generally good and don't show major discrepancies. The 
likely explanation for this difference is that the model ionic 
abundances report ed here are fo r the entire nebula, while 
those reported by iTsamis et al.l (|2003r ) only apply to the 
parts of the nebula caught under the slit. 

Table [^ provides complete fractional ionic abundances 
for our NGC 6302 model for both the circumstellar disk and 
bipolar lobes. Figure [S] shows the ionic fractions for a num- 
ber of elements along two sightlines that pass through the 
bipolar lobes (along the z-axis) and the circumstellar disk 
(along the x and y axes). In the bipolar lobes nearly all ion- 
ization stages are seen somewhere, with a clear ionization 
sequence from the inner heavily ionized regions to the outer 
parts of the lobes where some material is neutral. Hydrogen 
is 73% ionized over the entire lobes, while helium is 51% 
singly ionized and 35% doubly ionized. In the circumstellar 
disk the ionization front occupies a thin region at the inner 
edge of the disk, almost one hundred times closer to the cen- 
tral star than in the bipolar lobes. Because of this the vast 
majority of the circumstellar disk is neutral with hydrogen 
>99% ne utral in the disk. This i s in agreement with the Hq 
images of iMatsuura et al.l (|2005r ) that show very little emis- 
sion from the disk of NGC 6302. Compari son between the 
[Ne v] images presented by iBohigaj (IQOS) and our model 
Ne*""" fractions also shows good agreement, with 13% of the 
neon in the bipolar lobes in the form of Ne*^ , but less than 
0.1% in the circumstellar disk. 

Beyond the H^ ionization front at a depth of '^ 1.4 x 
10^^ cm the disk is optically thick and our model had diffi- 
culties converging due to the lack of photons. This is evident 
from the noise in the ionic fractions, predominantly of neu- 
tral species, beyond the ionization front. However, we found 
very little variation between Monte Carlo simulations in the 
total neutral and singly-ionized fractions of the majority of 
species, indicating that this noise did not affect our predicted 
emission line fluxes. 



5.4 The nebular temperature structure 

Mean electron temperatures weighted by ionic abundance 
have been calculated and are given in Table 1101 The ionic 
temperatures for the same ions in different regions of the 
model show the influence of the density on the shielding 
of ions and therefore the temperature structure. In the op- 
tically thick circumstellar disk the neutral species all have 
temperatures 1000-3000 K, while those species in the op- 
tically thin bipolar lobes have temperatures around 7,000- 
10,000 K. The ionic temperatures increase towards higher 
ionization stages, but increase faster in the disk than in the 
lobes. 

Because mocassin does not currently treat the neutral 
region fully, the temperatures of the neutral species likely 
represents the narrow transition region between ionized and 
neutral regions and not the highly optically thick and neu- 
tral region deep in the circumstellar disk. Beyond the ion- 
ization front, temperatures in the circumstellar disk drop 



Table 11. Predicted unreddened, reddened, and observed pho- 
tometry of the central star of NGC 6302. All magnitudes are 
presented on the Vega system . Reddened magnitudes were calcu- 
lated using a iHowarthI 1 119831) extinction law with ij = 3.1 and 
Av = 8.22. 



Filter 


A 

(A) 


FWHM 

(A) 


Observed 


Model 




mo m 


F469N 
F673N 


4688.1 
6765.9 


37.2 
100.5 


22.63 ±0.10 
20.08 ± 0.07 


12.57 22.60 
13.73 20.05 



to well below 1000 K and some cells have temperatures 
< 100 K. In these co nditions comp l ex mo lecules are likely to 
form, as observed bv lPeretto et al.l (120071 ). but we are unable 
to model this with the current version of MOCASSIN. 

Electron temperature profiles are shown in Figure]^ as 
a function of radius along two axes. In the inner parts of 
the lobes the temperature reaches in excess of 5 x 10* K 
(due to the lack of neutral H'' which would otherwise cause 
Ly-Qf cooling) but drops to 10,000 K by a radial distance of 
~ 7 X 10^^ cm. The temperatures in the disk do not reach 
such high levels, but do show the characteristic drops at 
the He^"'", He""", and H"*" ionization fronts (the latter pair 
of which, at a radius of --^ 1.4 x 10^^ cm, have coincident 
ionization fronts as would be expected for such a hot central 
star) . This range of temperatures is in full agreement with 
that found by many previous authors (e.g. TableO, namely 
a nebula with a mean ion- weighted electron temperature of 
19,400 K in the bipolar lobes, 3,800 K in the circumstellar 
disk and 12,000 K averaged over the entire nebula. 



5.5 Comparison with observations of the central 
star 

The central star of NGC 6302 was first detected in HST 
observations by ISzvszka et al.l ([20091, and Figure [T|, who 
attempted photometry in a number of narrow filters. Due 
to nebular contamination in many of the filters only two 
of these provide a reliable measurement of the stellar mag- 
nitude. Table [TT] provides a comparison of their measured 
values with predictions from our model. To calculate the 
unreddened stellar magnitudes the T^ff = 220, 000 K, 
Li, = 14, 300I/Q model atmosphere was convolved with the 
relevant HST filter profiles and a spectrum of Vega. The 
resulting magnitudes are ~6-10 magnitudes brighter than 
the observed phot ometry suggesting a c onsiderable amount 
of extinction (e.g. IMatsuura et al.ll2005h. Us i ng an R = 3.1 
extinction law in the form given bv iHowarthI (|l983r ) we red- 
dened the model photometry, varying the extinction to ob- 
tain agreement with the observed magnitudes. A very good 
fit was found using Ay = 8.22, fitting both magnitudes 
to within the observational uncertainties. T his extinction 
is slig htly larger than the Ay ~ 6.7 found bv ISzvszka et al.l 
( 20091), but in a greement with other literature values (e.g. 
IMatsuura et al.i2005 . estimated Ah a = 5 — 7, approximately 
Av =6—8). The good agreement between observed and pre- 
dicted photometric magnitudes and colours is strong verifi- 
cation of the central star luminosity and effective tempera- 
ture we derived. 

As illustrated by Figure [5] the high temperature and 
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Table 9. Nebular-averaged fractional ionic abundances for the NGC 6302 model. For each element the upper row is for the optically 
thick circumstellar disk and the lower row is for the optically thin bipolar lobes. 













Ion 














Element 


I 


II 


III 


IV 




V 


VI 




VII 




VII 




H 


0.991 
0.272 


0.937(-2) 
0.728 






















He 


0.996 

0.142 


0.369(-2) 
0.506 


0.236(-3) 
0.352 




















C 


0.983 


0.150(-1) 


0.157(-2) 


0.821 


-4) 


0.408(-4) 


0.255 


-6) 


0.504 


-9) 








0.366(-2) 


0.497 


0.255 


0.910 


-1) 


0.151 


0.281 


-2) 


0.346 


-4) 






N 


0.982 


0.179(-1) 


0.994(-4) 


0.502 


-4) 


0.543(-4) 


0.243 


-4) 


0.177 


-7) 


0.002 


-9) 




0.137 


0.479 


0.135 


0.495 


-1) 


0.710(-1) 


0.129 




0.280 


-3) 


0.269 


-6) 


O 


0.981 


0.187(-1) 


0.111(-3) 


0.501 


-4) 


0.588(-4) 


0.199 


-4) 


0.235 


-5) 


0.159 


-9) 




0.311 


0.316 


0.136 


0.391 


-1) 


0.737(-l) 


0.935 


-1) 


0.317 


-1) 


0.468 


-5) 


Ne 


0.983 


0.153(-1) 


0.105(-2) 


0.443 


-4) 


0.803(-4) 


0.996 


-5) 


0.260 


-6) 


0.229 


-7) 




0.196(-1) 


0.210 


0.528 


0.365 


-1) 


0.134 


0.671 


-1) 


0.422 


-2) 


0.851 


-3) 


Na 


0.883 


0.116 


0.133(-2) 


0.583 


-4) 


0.727(-4) 


0.196 


-4) 


0.399 


-5) 


0.732 


-6) 




0.315(-1) 


0.260 


0.450 


0.473 


-1) 


0.878(-l) 


0.681 


-1) 


0.411 


-1) 


0.114 


-1) 


Mg 


0.736 


0.264 


0.284(-3) 


0.213 


-4) 


0.663(-4) 


0.410 


-4) 


0.753 


-5) 


0.119 


-5) 




0.271(-1) 


0.586 


0.140 


0.205 


-1) 


0.635(-l) 


0.904 


-1) 


0.529 


-1) 


0.167 


-1) 


Al 


0.247 


0.752 


0.512(-3) 


0.190 


-4) 


0.653(-4) 


0.435 


-4) 


0.107 


-4) 


0.171 


-5) 




0.162(-2) 


0.612 


0.123 


0.334 


-1) 


0.601(-1) 


0.806 


-1) 


0.623 


-1) 


0.220 


-1) 


Si 


0.249 


0.750 


0.403(-3) 


0.421 


-4) 


0.555(-4) 


0.195 


-4) 


0.653 


-5) 


0.172 


-5) 




0.134(-2) 


0.734 


0.215(-1) 


0.458 


-1) 


0.510(-1) 


0.600 


-1) 


0.561 


-1) 


0.248 


-1) 


S 


0.690 


0.305 


0.472(-2) 


0.140 


-4) 


0.427(-4) 


0.561 


-4) 


0.237 


-4) 


0.364 


-5) 




0.597(-3) 


0.390 


0.341 


0.324 


-1) 


0.331(-1) 


0.609 


-1) 


0.966 


-1) 


0.373 


-1) 


CI 


0.911 


0.862(-l) 


0.255(-2) 


0.722 


-4) 


0.395(-4) 


0.292 


-4) 


0.195 


-4) 


0.432 


-5) 




0.221(-2) 


0.447 


0.244 


0.931 


-1) 


0.262(-l) 


0.359 


-1) 


0.857 


-1) 


0.558 


-1) 


Ar 


0.967 


0.325(-l) 


0.722(-3) 


0.344 


-4) 


0.101(-4) 


0.514 


-4) 


0.687 


-4) 


0.591 


-5) 




0.309(-l) 


0.198 


0.432 


0.973 


-1) 


0.121(-1) 


0.405 


-1) 


0.135 




0.464 


-1) 


K 


0.874 


0.125 


0.654(-3) 


0.618 


-4) 


0.271(-4) 


0.402 


-4) 


0.419 


-4) 


0.782 


-5) 




0.482(-l) 


0.387 


0.265 


0.750 


-1) 


0.173(-1) 


0.309 


-1) 


0.920 


-1) 


0.683 


-1) 



luminosity of the central star combined with its hydrogen- 
deficient nature results in a not-insignificant flux at UV 
and X-ray energies such that the central star of NGC 6302 
should be a detectable X-ray source. This provides an ob- 
servable test of both the high temperature and hydrogen 
deficiency of the central star. We have calculated the lumi- 
nosity of the central star in the Chandra (0.5-8.0 keV) and 
XMM-Newton (0.3-4.5 keV) bands as 5 x 10^" erg s"^ and 
6 X 10''^ erg s~^, respectively. X-ray absorption from neu- 
tral hydrogen in the circumstellar disk is calculated using 
the extinction of Av = 8.22 estimated above , and the con - 
version to X-ray absorbing cross-section from lRvterl (|l996r ). 
At a distance of 1.17 kpc this predicts absorbed fluxes of 
~ 5 X 10^^^ erg s^^ cm^^ for both X-ray bands (absorption 
due to neutral hydrogen is stronger at lower energies, off- 
setting the lower energy range of the XMM-Newton band 
for this soft X-ray source). Unfortunately the uncertainty 
on the predicted flux is relatively high due to uncertainties 
on the exact form of the stellar spectrum and the column 
of absorbing gas along the line of sight. However a H-rich 
central star of an equivalent effective temperature would be 
a signiflcantly weaker X-ray source such that the detection 
of an X-ray point source in NGC 6302 would be sufficient to 
confirm the H-deficient nature of the central star. Since X- 
rays may also originate from shocked ga s or collimated flows 
in the central cavi ty of the nebula (e.g. iKastner et al.ll2008l : 
IParkin et al.ll2009l ). such observations must also distinguish 



between diffuse and point-source X-ray emission if the X-ray 
detection of the central star is to be observationally tested. 



6 DISCUSSION 

In this section we discuss the properties of the nebula and 
central star modeled here with particular focus on the parts 
of the model that are most strongly constrained by the ob- 
servations. We discuss these in the context of previous stud- 
ies of NGC 6302 and the general properties of planetary 
nebulae. Finally we suggest an evolutionary scenario for the 
formation and history of NGC 6302. 



6.1 The nebula structure 

NGC 6302 has long been known as a highly inhomogeneous 
nebula and our model matches this. Between the bipolar 
lobes and the inner edge of the circumstellar disk we model a 
density contrast of two orders of magnitude, from 2000 cm~'^ 
in the lobes to 300,000 cm~^ in the circumstellar disk. This 
was necessary to reproduce many of the diagnostic emission 
line ratios, which were impossible to fit without including 
both the low and high density par ts of the model. Such 
high densities are rare in PNe (e.g. lOsterbrock fc Ferlandl 
liooe). In their study of FN densities 'Tsamis ct al." (200^, 
NGC 6302 consistently displayed the highest densities from 
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Figure 6. Ionic fractions and electron temperatures for two sightlines tiirough the model nebula. Left: along the z-axis through the 
bipolar lobes; right: along the x- or j/-axis through the circumstellar disk. Note the different spatial dimensions on each. Ionization 
fractions are shown for helium and hydrogen (2nd panel from top), carbon (middle panel), oxygen (2nd panel from bottom), and neon 
(bottom panel). The dashed line in the right-hand panels indicates the edge of the circumstellar disk. 



Table 12. Model properties of the best fitting final model of 
NGC 6302. 



Model property 




Value 


Uncertainty 


Volume of bipolar lobes 
Mass of bipolar lobes 
Volume of circumstellar disk 
Mass of circumstellar disk 


9.30 X 
1.40 X 


10^3 cm^ 

2.47 M© 
10^3 cm3 
2.22 Mq 


±20% 
±30% 
±50% 
±50% 


Total modeled nebular mass 
Total ionized mass 




4.69 M© 
1.82 M© 


±40% 
±20% 



a range of diagnostic line ratios, with only NGC 5315 and 
IC 4191 showing comparably high densities from their sam- 
ple. 

Derived total volumes and masses for our model nebula 
are listed in Table [T^l The total ionized mass of our model 
is ~1.82 M© (predominantly located in the bipolar lobes), 
much higher than the typical ionized masses of Typ e I PNe 
(e.g. 0.54 M© for six MageUanic Cloud Type I PNe. iBarlowl 
[l987|). The total mass of nebular material modeled (not in- 
cluding dust) is ~4.7 M©, also notably higher than the typi- 



cal masses of PNe (jPottaschll 19961 . studied 96 PNe and found 
masses in the range 0.01 < Mi/M© < 0.8), but not implau- 
sible. The final disk radius value of 3 x 10^^ cm was con- 
strained from fitting the infrared SED (not discussed here), 
but is unconstrained from photoionization modeling. The 
total disk mass is in good agreement with previous observa- 
tions. IPeretto et al.l ()2007l ) estimate the mass of the molec- 
ular region of NGC 6302 to be as much as 1.8 ± 0.9 M©, in 
good agreement with our disk mass o f ^2.2 M©, of which 
97% (2.1 Mq) is neutral. However, IPinh-V-Trung et all 
(J2008I ) derive a significantly lower molecular mass of 0.12 M© 
(scaled to a distance of 1.17 kpc), but note that this may 
be a lower limit if the CO temperature in the disk is higher 
than the 30 K they assume. 

During the modeling process we attempted to model the 
circumstellar disk as a torus but found that the height of the 
circumstellar material was well constrained by the observa- 
tions and that a fiattened circumstellar disk was a better fit. 
Evidenc e for a n extended disk- like morphology comes from 
iBohigasI ([199J) who observed differential extinction in the 
western lobe, whi ch he sugges ted is viewed through the outer 
parts of the disk. Icke (2003) suggested that a warped disk 
structure is also necessary to form the multiple lobes seen in 
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Table 10. Mean electron temperatures (K) weighted by ionic species for the NGC 6302 model. For each element the upper row is for 
the optically thick circumstellar disk and the lower row is for the optically thin bipolar lobes. 











Ion 








Element 


I 


II 


III 


IV 


V 


VI 


VII 


VIII 


H 


2018 
7742 


7367 
12682 














He 


2237 
7468 


8094 

8874 


14838 
15415 












C 


896 


4736 


5754 


15740 


16164 


17461 


20243 






8088 


8164 


9769 


15372 


18542 


22326 


36914 




N 


1177 


4817 


11435 


15701 


15813 


16424 


18225 


20780 




7259 


8478 


10270 


13823 


16863 


18911 


24343 


43851 


O 


1428 


4596 


10912 


15707 


15807 


16330 


18343 


20485 




7681 


8764 


10358 


14707 


16787 


18489 


20652 


29473 


Ne 


944 


4421 


8724 


15628 


15861 


16720 


19346 


21347 




6964 


7740 


9147 


14256 


17270 


19222 


22298 


28040 


Na 


2041 


4061 


8681 


15474 


15785 


16139 


17509 


19813 




7804 


7818 


9164 


13299 


16492 


18159 


19417 


21362 


Mg 


2002 


3720 


7517 


15648 


15750 


15945 


16777 


18939 




8602 


8358 


9933 


13272 


15950 


17627 


19022 


20724 


Al 


2785 


2621 


6110 


15370 


15742 


15893 


16545 


18455 




8533 


8452 


9680 


11991 


15757 


17382 


18807 


20286 


Si 


2806 


2724 


8260 


15705 


15753 


16116 


17115 


18994 




8720 


8718 


11983 


13996 


16281 


17937 


18963 


20207 


S 


440 


3389 


4853 


15398 


15715 


15781 


16219 


17532 




7954 


7967 


9250 


11327 


14722 


16500 


18216 


19562 


CI 


1112 


4072 


4904 


14724 


15693 


15801 


16278 


18064 




7986 


8057 


9390 


11455 


15087 


16620 


18059 


19445 


Ar 


848 


4356 


8356 


10308 


15681 


15732 


15911 


16941 




7037 


7624 


8789 


10223 


13535 


15601 


17615 


19354 


K 


2113 


4509 


9048 


13819 


15641 


15738 


15946 


17080 




8079 


8046 


9385 


10977 


14358 


15994 


17553 


19050 



NGC 6302. A thin gas disk and a molecular torus could both 
be present in NGC 6302 if the torus of dust and molecular 
material were embedded within the outer regions of a suffi- 
ciently flared circumstellar disk. It is possible that any once 
torus-like structure would become thin and disk-like due to 
radiation pressure and stellar winds. 



6.2 Nebular abundances 

We fit our model of NGC 6302 with a helium abun- 
dance of 0.153, which is high even for a Type I PN 
IjKingsburgh fc Barlowlll994 ). indicative of a large amount 
of stellar processing in the progenitor star. According to 
the dredge-up models of lBecker fc IberJ (| 19801 ). nebular ma- 
terial with this He abundance should correspond to a star 
with an initial mass > 4 M© . The O/H abundance we derive 
for NGC 6302 (5.1 ± 1.3 x 10"'') agrees within the uncer- 
tai nty limits with the mean O/H ratio for Galactic disk PNe 
bv lKingsburgh fc Barlow ll 19941) and with the solar O/H ra- 
tios of lCaffau et al.l (|200a ) and lAsplund et all (|2009l ). We do 
not confirm the deplet ion of oxygen by a fact or of two that 
was fo und bv lPottasch fc Beintemal (|l999r ) and lTsamis et al.l 
(|2004l ) from empirical abundance analyses of coUisionally 
excited UV, optical, and infrared lines. We attribute this 



difference to the neglect by empirical analyses of the large 
fraction of oxygen that resides in very high ionization stages 
in this unusually highly ionized nebula (see Table [9]). 



We find a C/0 ratio of 0.43 for NGC 6302, indicat- 
ing that the nebula is predominantly 0-rich, in agreement 
with all other measuremen ts of the C/O ra tio in the lit- 
erature (e.g. C/0= 0.20. lAller et al.l ll98]J. C/0= 0.26 
IPott asch fc Beintemal Il999l. C/0= 0.88, ICasassus et^ 
200d. and C/0= 0.31. iTsamiset "all l2003l ). with the ma- 
jority of values < 0.5. Our derive d C/H ratio of 2.2 x lO""' 
is 80% of the solar value listed bv lAsplund et al.l l|2009l ). in- 
dicating only a slight depletion of carbon, if any. However, 
the N/H ratio we derive for NGC 6302 (3.9 x 10"'') is a fac- 
tor of 5.7 larger than the solar value given bv lAsplund et al.l 
(120091) ■ too large to have been produced by secondary con- 
version of initial carbon, while the combined (C-|-N-|-0)/H 
ratio for NGC 6302 of 11.2 x 10"'' is 35% larger than the solar 
value of 8.28 x 10"*. We conclude that primary enrichment 
of nitrogen occurred in the precursor star of NGC 6302, via 
3rd dredge-up enhancement of carbon, followed by hot bot- 
tom burning CN-cycle conversion of dredged-up carbon to 
nitrogen. 



Several elements appear particularly depleted, such as 
magnesium (by a factor of two compared to the solar 
value) and aluminium (by a factor of one hundred). This 
is strong evidence for depletion onto dust grains, and given 
the evidence for la rge amounts of dust in NGC 6302 (e.g. 
iMolster et al.ll200ll ). is not unexpected. 
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6.3 The central star 

The observation of emission lines from very high ionization 
stages and the nebula's Type I status has long suggested 
that the central star of NGC 6302 was extremely hot. The 
hottest PN central stars directly measu red to date have sur - 
face temperature of up to 180,000 K (jWerner et al.lll997f) . 
which though high, are lower than all e stimated temper- 
ature s for NGC 6302's central star (e.g. lAshlev fc Hvlandl 
119881 . estimated a temperature of 430,000 K). 

Using a H-deficient stellar model atmosphere we find a 
best fit to the observed emission-line spectrum with cen- 
tral star properties Tefj — 220,000 K, \ogg — 7.0, and 
L* — 14,300 L0. This is the lowest temperature ever esti- 
mated for the central star of NGC 6302, due to our use of a 
H-deficient stellar atmosphere which has a larger flux at the 
energies necessary to significantly populate the highest ion- 
ization states from which emission lines are observed (see 
Figure [5]). However, this luminosity is significantly larger 
than previous estimates, or our estimate of the total lumi- 
nosity of 5700 Lq for NGC 6302, in Section [231 from in- 
tegrating its dereddened observed emission across all wave- 
lengths longwards of 1200 A. That this is only 40% of the 
stellar luminosity needed by our bipolar model to match the 
observed line and continuum fluxes suggests that 60% of the 
stellar luminosity escapes through optically thin regions of 
the nebula. 

The observed expansion of NGC 630 2's bipolar lobes 
JMeaburn et al.l l2008l : ISzvszka et al.l l201ll 'l suggests it once 
had a fast wind that could affect the stellar spectrum. How- 
ever strong mass-loss is usually only seen at lower effec- 
tive temperat ures and lower sur face gravities, e.g. log;? = 
3.0 — 5.5 (e.g. iBohlin et al.lll982l . find no evidence for stel- 
lar winds in the UV spe ctra of very hot PN central stars 
that [Benedict et al.ll200 9 found to have surface gravities of 
log gr = 6.9 — 7.3). As we will show in Section [6.41 evolution- 
ary tracks for LTP and VLTP models imply a surface gravity 
of log (7 ~ 6.6 for the Te// and L* we have derived, and this 
can also be derived from simple stellar relat ions by assum- 
ing a typical central star mass. Furthermore ICasassus et al.l 
(|2000f ) could find no evidence for a hot wind or a wind blown 
cavity from observations of the IR coronal lines originating 
in the center of NGC 6302. 

If the central star of NGC 6302 is H-deficient this im- 
plies that it has undergone some sort of late thermal pulse 
during its evolution, which has co nsumed the majorit y of the 
remaining surface hydrogen (e.g. lHerwig et al.lll999l ). There 
are three possible scenarios for this: the 'after final thermal 
pulse' (AFTP) scenario, where an e xtra thermal p ulse is ex- 
perienced at the top of the AGB (|Herwiell200ll ): the 'late 
thermal pulse' (LTP) that occurs during the post- AGB evo- 
lution while H-burning is still active (Blocker 2001); or the 
'very late thermal pulse' (VLTP) exp erienced by a hot whit e 
dwarf during its early cooling phase (|lben fc Renzinilll983h . 
Only the VLTP scenario consumes the majority of the re- 
maining hydrogen in the central star, while the other two 
scenarios suggest a dilution of the surface hydrogen to very 
low levels. We are currently unable to compare model at- 
mospheres with H-deficient abundances with H-poor abun- 
dances, so we cannot use this method to separate these sce- 
narios. 

H-deficient central stars include the Wolf-Rayet central 



stars ([WC], e.g. Hamann 19971) and th e hotter and rarer 
PC 1159 stars l| Werner fc Herwid 1200 6*). [WC]-type stars 
have been sugg ested to be the evolutio nary precursors of 
PC 1159 stars l|G6rnv fc Tvlendal I2OO0I ). with the earlier 
stages characterised by strong stellar winds and interaction 
with the remaining nebulous material, and the latter stage 
exhibiting an often isolated central star at its hottest tem- 
perature. The high temperature and high surface gravity 
of the central star of NGC 6302 suggests it is similar to a 
PC 1159 star, and the presence of the remaining nebulous 
material can be explained if the star were particularly mas- 
sive and has therefore evolved to this temperature before 
the nebula has dispersed. 

There are many similarities between NGC 6302 and 
other PNe with H - defic ient central stars that support this 
assertion. Zijlstra, ( (200 1 |) notes that nearly all PNe with 
[WC]-type central stars feature both O-rich and C-rich 
circumstell ar material (i.e. dual-dust chem istry), as does 
NGC 6302 (|Molster et al.l200ll : [Cohenl200ll ). However, most 
H-deficient nebulae exhibit much stronger PAH (polycyclic 
aromatic hydrocarbon) emission than NGC 6302, though 
this could be due to the very low C/O ratio in NGC 6302 
(e.g. iGuzman-Ramirez et al.l I2OIII . show that the produc- 
tion of PAHs is greatly reduced at low C/O ratios). An- 
other similarity between NGC 6302 and PNe with H-poor 
or H-deficient central stars is that they both show highly 
inhomogeneous s tructures with many kn ots and filaments 
(see Figure [Tl and iGornv fc Tvlendall200d 'l. Such knots may 
be associated with th e complex turbulent velocities found 
bv I Acker et~ai] (|2002l ') in the lobes of the majority of [WC]- 
type PNe. However the k nots observed in NG C 6302 do not 
appear to be H-deficient (jSzvszka et al.ll201lh . unlike those 
observed in other H-poor PNe (e.g. Abell 30 and Abell 78). 
This could indicate that the knots in NGC 6302 were ejected 
during an earlier, H-rich phase of mass-loss and have since 
been disrupted and accelerated by a later, H-poor mass-loss 
event. 



6.4 Evolutionary tracks for late thermal pulse 
evolution 

To further study the central star of NGC 6302 we have com- 
pared its properties determined from our models with those 
from evolutionary tracks for late thermal pulse models. This 
might allow us to determine the core mass of the central 
star, and therefore potentially its original mass, as well as 
its evolutionary state. 

Figur e [7| shows evo lutionary tracks for the VLTP 
model s of iBldckerl (119951) and the LTP and VLTP mod- 
els rfjMmer^ertoiamr^Aithaui (|2006l ). In the models 
of iBlockeiJ (|l995l ) the central star of NGC 6302 falls 
between the 0.61 and 0.84 M© core mass tracks, at 
^0. 82 M(7), corresponding to an initial mass of 4.8 Mq. In 
the iMiller Bertolami fc Althaud (|2006l ) models the central 
star corresponds to a star with a core mass of ~0.73 M©, an 
initial mass of 3.7 Mq, and log g = 6.6, close to our modeled 
value. 

We can compare the timescales of these different 
evolutionary tracks with the known evolutionary history 
of NGC 6302, such as the last major mass-loss event 



2200 yrs ago (JMeaburn et al.1 I2OO8I : ISzvszka et al] I2OIII ). 
The timescales for LTP events are considerably longer than 
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Figure 7. Hertzprung-Russell diagrams for very late thermal pulse evolutionary models co mpared to the position of the central star of 
NGC 6302 derived from our models (with 20% error bars). Left: Evolutionary tracks from iBlocken ||l995|) with core masses indicated. 
Tracks show the first PN central st ar phase, the born-again event and the second PN central star phase. Right: Evolutionary tracks from 
iMiller Bertolami &: Althaug Il200y) with core masses (indicated) showing only the second PN central star phase. 



those for VLTP events, which could suggest that the VLTP 
event is more likely given the short timescale since the last 
mass- loss event. Consider i ng th e VLTP model timescales, 
for the models of iBlockerl (119951 ) for a 4.8 Mq initial mass 
star the position corresponds to an age of ~2100 yrs since 
the star left the end of the AGB and ~300 years since the 
VLTP event, suggesti ng that the eruptive event traced by 
iMeaburn et al.l (|2008h was the last mass-loss event on the 
AGB before the ce ntral star evolved to hotter tempera tures. 
For the models of iMiller Bertolami fc AhhausI l|2006l ) only 
ages since the VLTP event are provided, which, for a 3.7 Mq 
initial mass star correspond to ~2000 yrs since the helium 
flash that caused this event, consider ably longer th an the 
equivalent evolutionary period in the iBlockerl (|l995l ) mod- 
els. 

If the timescales of IMiller Bertolami fc Althauj (120061 ) 
are correct this s ugges ts that the eruptive event identified by 
iMeaburn et al.l (|2008l ) is associated with a VLTP event. A 
VLTP mass - loss event however is not thought to be partic- 
ularly large. Ivan Hoof et al.l (|2007l ) estimated that the mass 
lost in the VLTP eruption of V334 Sgr might be as much 
as 0.01 Mq, at least an order of magnitude smaller than 
would be necessary for one pai r of NGC 6302's bipolar lobes. 
However. lAlthauset all (|2008l ) found evidence that PC 1159 
stars may have much thinner He-rich envelopes than previ- 
ously thought and this might be due to considerable mass- 
loss during the VLTP event, though the helium layer is only 
~0.02 Mq thick and this pu ts an upper limit on what can 
be lost during such an event. ISzvszka et al.l l|201H ) note that 
the mass and momentum of the bipolar lobes is too high for 
them to have been ejected in a single mass-loss event and 
suggest that the recent mass-loss event may have accelerated 
gas that was previously ejected. 

The large discrepancy in timescales between the models 
presented here for similar VLTP events implies a significant 
level of uncertainty in the current models that requires more 
accurate evolutionary tracks. Therefore at this point it is 



impossible to distinguish between the two sets of models 
as either could imply a timescale that agrees well with the 
evolutionary history of NGC 6302. Combining the derived 
core mass (0.74-0.82 Mq) with the nebular mass estimated 
from this work (4.7 Mq; Table [12)) implies an initial stellar 
mass of at least 5.5 Mq. The initial stellar mass of 3.7- 
4.8 Mq from evolutionary tracks is not far from the 5.5 Mq 
that we derive. 



6.4.1 Comparison with other known PG 1159 stars 

A number of PG 11 59 stars have measured e ffective tempera- 
tures and gravities (| Werner fc Herwig|l2006l ). Figure [S] shows 
these stars in the log g - log Te// plane (chosen because these 
quantities are better constrained than their luminosities), 
alongside the central star of NGC 6302 (using the value of 
logg ~ 6.6 fitted from Figure [T] which is better constrained 
than our modeled value), an d the helium-deficient star 
H1504-I-65 (|Werner et al■ll2004^■ In addition we show evo - 
lutionary tracks for the IMiller Bertolami fc AlthausI (|2006l ) 
LTP/VLTP models. The position of the central star of 
NGC 6302 is in good agreement with that of many of the 
other PG 1159 stars. Is is clearly hotter than any of the 
known PG 1159 stars, and even hotter than the previously 
hottest known central star, H1504-f 65, though that source 
does appear to have had a higher initial mass and a previ- 
ously higher temperature (it is now on the cooling track). 



6.5 The evolutionary history of NGC 6302 

From the evidence presented here we suggest that NGC 6302 
was created by a star with an initial mass of at least 
5 Mq, in agreement with both the chemical abundances 
and evolutionary tracks for the central star. The star 
may have had a lower-mass binary companion that di- 
verted the AGB-phase mass-loss from isotropic to equato- 
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Figure 8. VLTP evolutionary trades showing the second PN cen- 
tral star phase in t he logq - logT^ff plane with core masses in- 
dicated llMiller Bertolami fc Althauj boOSh . The central star of 
NGC 6302 is shown with 20% error bars. The positions of the 
known PNe with PG 1159 ty pe central stars are show n as open 
circles with data taken from IWerner fc Herwig 11200611 and the 
helium-deficient star H1504+65 llWerner et alj l2004h is shown 
with a cross. 



The circumstellar disk would h ave caused later mass- 
loss events to form the bipolar lobes. iMatsuura et al.l ([20051) 
observed the circumstellar disk to be warped, which may 
be caused by an interaction with a bin ary companion, and 
according to the model of llckd (|2003l ) would cause each 
pair of bipolar lobes to form with a different axis, as is ob- 
served. Hot-bottom burning at the base of the convective 
envelope (which occurs for stars more massiv e than ~4 M0 , 
llben fc Renzinill 19831 : iBoothrovd et al.ll 19931 ) . would convert 
carbon to nitrogen and prevent the star or the ejecta going 
through a prolonged C-rich phase. 

At some point the star experienced some sort of late 
thermal pulse, ejecting a very small amount of material, 
causing the central star to become H-deficient and returning 
to the tip of the AGB. It is not currently possible to ob- 
servationally distinguish between the different LTP scenar- 
ios, t hough the evolutionary timescales for the VLTP mod- 
els of lMiller Bertolami fc Althaui (|2006l ') show the strongest 
agreement with the known mass-loss history. Both sets of 
evolutionary tracks (above) suggest that the temperature of 
the star will continue to rise, potentially reaching a maxi- 
mum temperature of ~280,000 K, before entering the white 
dwarf cooling track again. 



LividI 19901 : 



rial, forming the circumstellar dis k (e.g. lBond 
iGarcia-Arredondo fc Frankll2004l 'l. 

The question of the age and stability of the circumstellar 
material is important for understanding t he history of such a 
nebula. The presence of large dust grains (IHoare et al.lll992l) 
and a high fraction of crystalline material i Molster et al.l 
[2OOII) could indicate that the circumstellar disk is long- 
lived, providing a shielded and stable environment for these 
grains to build up. However, such dust characteristics could 
also be caused by high mass-loss rates and densities in the 
ejected material (ISvlvester et al.l Il999l . found that the de- 
gree of crystalli nity in OH/IR stars appears to scale with 
mass- loss ratel . IPeretto et al.l ((2007]) measure an expansion 
velocity of ~8 km s~^ for the circumstellar torus and argue 
that it was ejected during a ~5000 year mass-loss event that 
started ~7500 years ago. Such a scenario requires mass-loss 
rates of 5 x 10~* M© yr~^ over this time period, potentially 
high enough to create the conditions for grain growth and 
crystallization in NGC 6302, but higher than current evolu- 
tionary models can sustain, even for shorter periods of time 
(e.g. Vassihadis fc Wood 1993; Blocker 1995). 

The kinematic age may be a lower limit if the disk has 
been accelerated by radiation pressure from the highly lu- 
minous central star. A simple calculation for the momentum 
imparted by radiation pressure on the circumstellar disk sug- 
gests that, based on the param eters derived here an d the 
expansion velocity measured by IPeretto et al.l ((20071), this 
would take ~500,000 years at the current stellar luminos- 
ity (assuming that ~10% of the stellar luminosity intercepts 
the circumstellar disk). This timescale is significantly longer 
than typical PNe lifetimes and as such unfeasible as the ma- 
jor acceleration mechanism fo r the circumstellar di sk. A sim- 
ilar argument is presented bv IPeretto et al.l ((2007]) who also 
conclude that alternative mechanisms are required and sug- 
gest interaction with a binary companion as a likely propul- 
sion mechanism. 



7 CONCLUSIONS 

We have constructed a 3D photoionization model of the ex- 
treme planetary nebula NGC 6302 using the Monte Carlo 
photoionization and radiative transfer code MOCASSIN. Our 
model reproduces the majority of emission line strengths 
and ratios, which place strong constraints on many of the 
stellar and nebular properties. The model consists of a very 
dense circumstellar disk where densities reach 300,000 cm~"^ 
at the inner edge, a large pair of bipolar lobes with a con- 
stant density of 2000 cm~^, and an intermediate component 
we have dubbed 'the outflow'. This combination of compo- 
nents was required to match the majority of density- and 
temperature-sensitive line ratios from which a wide range of 
densities and temperatures have been observationally deter- 
mined by previous authors. A number of line ratios remain 
not matched, which we suggest is due to complex density 
inhomogeneities throughout the nebula. 

We derive a total nebular mass of 4.7 M©, of which 
1.8 Mq is ionized, almost entirely in the bipolar lobes. The 
C/O ratio for NGC 6302 is 0.43 indicating a predominantly 
0-rich nebula, in agreement with all other measurements 
in the literature. Carbon is found to be marginally under- 
abundant compared to the solar value, while nitrogen is sig- 
nificantly over-abundant. The combined (C+N+0)/H abun- 
dance is 35% larger than the solar value, all of which suggests 
that a 3rd dredge- up occurred in the precursor of NGC 6302, 
enriching the central star in carbon and nitrogen, followed 
by hot bottom burning CN-cycle conversion of carbon to 
nitrogen. 

In modeling the central star we have incorporated 
NLTE model atmospheres to reproduce the ionizing flux. 
Fits to the optical emission-line spectrum imply an ex- 
tremely high temperature for the central star, and we used 
the observed emission from high ionization-stage infrared 
'coronal' lines to further constrain the form of the ioniz- 
ing flux distribution. Using a solar abundance stellar at- 
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mosphere we were unable to fit all the observed line fluxes, 
overestimating many of the lower ionization stages while un- 
derestimating the higher stages. A substantially better fit 
was obtained using a hydrogen-deficient stellar atmosphere, 
implying that the central star of NGC 6302 is likely to be 
H-deficient. 

Finally we compare the properties of the central star 
with evolutionary tracks for late thermal pulses that are ca- 
pable of removing the majority of hydrogen from the central 
star. We find a good fit to a very late thermal pulse track 
for a star with an initial mass of 4-5 M0. Timescales for 
this evolutionary model imply that the central star left the 
top of the AGB ~2100 years ago, in good agreement with 
the age of one of the pairs of bipolar lobes determined by 
iMeaburn et al.l (|2008l 'l to be ~2200 years. This mass is in 
reasonable agreement with the total modeled nebular mass 
plus central star mass of 5.5 M©, taking into account con- 
tributions from the dust component. The mass is also in 
agreement with the chemical abundances in the nebula such 
as the high helium abundance, low C/0 ratio and slightly 
enhanced (C-fN+0)/H ratio. 

A future paper will discuss the dust component of this 
model that used mocassin's ability to fully model the radia- 
tive transfer of dust in an ionized nebula. The model includes 
the majority of observed amorphous and crystalline dust 
species, using multiple dust chemistry distributions, non- 
standard grain size distributions, and a new light scattering 
code for non-spherical dust grains. 
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